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INTRODUCTION 


The  use  of  composite  materials  and  structures  to  provide 
characteristics  unattainable  directly  from  the  constituent 
materials  is  well  known.  Perhaps  the  most  widespread  example  is 
steel-reinforced  concrete  for  structural  applications  wherein  the 
high  tensile  strength  of  the  steel  in  conjunction  with  the  high 
compressive  strength  of  the  concrete  yields  a  composite  material 
with  structural  properties  far  superior  to  those  of  either 
component.  More  recently,  work  has  been  undertaken  to  apply  this 
principle  to  the  development  of  dielectric  materials  for  energy 
storage  applications.  It  is  anticipated  that  this  technology  will 
provide  materials  for  use  in  capacitors  with  greater  energy 
density,  lower  loss,  and  higher  breakdown  resistance.  This 
report  describes  the  results  of  an  effort  to  analytically  model 
the  net  or  resultant  dielectric  behavior  of  such  composite 
dielectric  materials. 

The  net  dielectric  behavior  of  a  randomly  interspersed 
composite  is  dependent  on  the  spatial  dimensionality  (1-D  vs.  2- 
D  vs.  3-D),  domain  geometries  (domain  size,  domain  shape, 
stratification,  etc.),  and  interconnection  effects  (percolation). 
This  collective  dependence  is  only  partly  understood.  One 
averaging  law  which  is  useful  to  the  experimentalist  because  of 
its  relative  ease  in  dealing  with  multiple  components  is  [1] 

ea  =  £  vk  e£  (1) 

where  e  is  the  net  dielectric  permittivity  of  the  composite  while 
is  the  dielectric  permittivity  component  occupying  volume 
fraction  Vj..  The  exponential  factor  a  depends  on  the  geometry  of 
the  constituent  components  and  has  been  rigorously  derived  for 
only  the  special  cases  of  layers  oriented  perpendicular  to  the 
applied  electric  field  and  layers  oriented  along  the  direction  of 
the  electric  field  [2].  For  the  case  of  layers  oriented 
perpendicular  to  the  applied  electric  field,  a  =  -1  and  equation 
1  may  be  interpreted  as  harmonic  averaging.  For  the  case  of 
layers  oriented  along  the  direction  of  the  applied  electric 
field,  a  =  +1  and  equation  1  may  be  interpreted  as  arithmetic 
averaging.  Although  no  obvious  physical  significance  may  be 
immediately  attached  to  the  case  of  a  =  0,  it  may  be  interpreted 
as  geometric  averaging.  No  generalized  analytic  expression  for  a 
exists,  and  a  primary  goal  of  this  work  is  to  develop  a  numerical 
technique  to  obtain  this  factor  for  different  composites. 

DIELECTRIC  PERMITTIVITY 

The  dielectric  permittivity  e  relates  an  applied  electric 
field  to  the  induced  displacement  field  arising  as  a  result  of 
the  applied  electric  field.  This  relationship  is  usually  written 


where  both  the  displacement  arid  electric  field  Ej  are  vector 


quantities  and  the  dielectric  permittivity  is  a  tensor  of 
rank  two  over  the  spatial  indices  i,j.  The  summation  convention 
applies  to  the  repeated  index  j.  In  the  case  of  isotropic  media, 
equation  2  takes  the  vector  form 

D  =  e  E  (3) 

with  e  now  written  as  a  scalar  quantity.  In  all  practical  cases, 
the  dielectric  permittivity  is  a  complex  quantity  and  may  be 
denoted 


e  =  e 


+  j  e" 


wherein  j=/-l.  The  Kramers-Kronig  relationship,  based  on  time- 
causality  considerations,  specifies  that  the  real  and  imaginary 
parts  of  e  are  not  independent.  In  engineering  applications,  the 
complex  nature  of  e  is  more  commonly  referred  to  by  the  use  of 
the  loss  tangent,  defined  as 


-  material  #1 


-  material  #2 


Figure  1.  Parallel  plate  capacitor  with  composite  dielectric. 


o 


Figure  3.  Four-terminal  network  pixel  element  model. 


Figure  2.  The  pixel  ~rid  forms  a  Cartesian  coordinate  system  with 
one  axis  parallel  and  one  axis  perpendicular  to  the  applied 
electric  field.  The  dielectric  properties  of  the  material  within 
any  pixel  are  assigned  by  means  of  the  pixel  grid  coordinate 
system. 

We  can  consider  each  pixel  site  as  a  four-terminal  network 
and  can  model  all  nearest-neighbor  interactions  of  any  pixel  by 
considering  the  array  of  nodes  formed  at  the  pixel  intersections. 
This  is  sufficient  to  assign  a  vector  displacement  to  each  pixel. 
The  four-terminal  network  within  each  pixel  is  comprised  of  a  set 
of  four  capacitors,  each  connected  between  a  common  node  at  the 
center  of  the  pixel  and  one  of  the  corners  as  illustrated  in 
Figure  3.  The  value  associated  with  each  capacitor  is  simply  the 
size-normalized  capacitance  between  the  pixel  center  and  any 
corner. 

We  can  redraw  the  composite  dielectric  pixel  grid  of  Figure 
2  with  the  four-terminal  equivalent  network  substituted  for  each 
pixel  element,  and  obtain  the  electrical  network  shown  in  Figure 
4.  This  network  may  be  analyzed  by  means  of  Kirchoff's  Laws  as 
applied  to  the  displacement  field  D  in  the  static  case  or 
displacement  current  dD/dt  in  the  dynamic  case. 

NETWORK  ANALYSIS 

In  the  application  of  Kirchoff's  Laws  to  the  analysis  of  the 
equivalent  electrical  network,  we  are  confronted  with  the  choice 
of  using  either  mesh  analysis  or  nodal  analysis.  We  choose  to 
employ  nodal  analysis  for  the  following  reasons: 

1)  The  numerical  solution  is  more  stable  using  nodal 
analysis  than  using  mesh  analysis.  Essentially,  as  one  traverses 
the  composite  network,  the  nodal  solution  transitions  smoothly 
from  pixel  to  pixel  and  the  relative  differences  in  currents  are 
small,  whereas  with  mesh  analysis  we  may  see  alternations  in  the 
sense  of  the  ] cop  currents  which  greatly  increases  the  relative 
differences  in  currents-  From  a  numerical  standpoint,  this 
choice  is  impcet-id'  in  redding  truncation  errors. 
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Figure  4.  Composite  dielectric  equivalent  electrical  circuit 
where  the  circuit  components  have  been  numbered  with  subscripts  1 
and  2  to  correspond  with  the  material  components  of  Figure  2. 
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2)  The  boundary  conditions  are  more  readily  implemented 
using  nodal  analysis.  This  includes  both  the  excitation  and 
lateral  boundaries. 

Implementation  of  the  nodal  analysis  results  in  a  succinct 
description  of  the  network.  The  currents  1^  between  the  ith 
node  and  its  nearest  neighbors  are  described  by 

2  Iij  =  0  (5) 

for  all  nodes  except  along  the  excitation  plane.  A  normalized 
current  is  applied  to  nodes  along  the  excitation  plane,  tor  which 
we  may  write 

2  Iij  =  j  w  (6) 

where  jw  is  the  complex  radian  frequency  of  the  excitation.  The 
summation  is  taken  over  the  nearest  neighbor  nodes.  The  size 
normalized  admittance  of  the  capacitor  elements  between  nodes  i 
and  j  may  be  written 

Yij  =  j  w  ei:j  (7) 

Equation  7  may  be  substituted  into  equations  5  and  6  using  Ohm's 
Law  and  the  potentials  at  nodes  i  and  j ,  resulting  in 

2  <4j  (Vi  -  Vj)  =  0  (8) 

for  most  nodes  and 

2  €ij  (Vi  -  Vj)  =  1  (9) 

for  nodes  on  the  excitation  plane.  The  solution  of  the  nodal 
equations  varies  with  frequency  since  the  of  the  constituent 

materials  are  in  reality  frequency  dependent  with  the 
different  materials  having  different  frequency  dependencies. 

BOUNDARY  CONDITIONS 

In  the  formulation  of  the  network  to  be  solved,  we  encounter 
two  common  types  of  Doundary  conditions  concerning  the 
displacement  field  at  the  lateral  boundaries.  The  first 
condition  is  that  of  'insulating'  sides,  wherein  we  consider  the 
dielectric  specimen  to  be  electrically  isolated  along  the  lateral 
boundaries.  In  this  case,  the  normal  component  of  the 
displacement  field  D  is  considered  to  be  zero  at  the  lateral 
boundaries.  For  simplicity,  we  also  assume  that  fringing 
electric  fields  are  nonexistent.  This  boundary  condition  is 
applicable  to  the  analysis  of  isolated  samples  or  lattice  cells 
with  mirror-symmetry-plane  boundaries.  The  second  condition  is 
that  of  'periodic'  boundaries,  wherein  we  consider  the  dielectric 
specimen  to  be  a  repeated  cell  in  a  cyclic  lattice  structure.  In 
this  case,  the  normal  component  of  displacement  field  is  non-zero 
and  the  same  at  both  boundaries,  as  is  the  tangential  component 
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Figure  5.  Lateral  boundary  conditions  on  the  composite 
dielectric. 


of  D,  The  two  types  of  boundary  conditions  are  illustrated  in 
Figure  5. 

COMPUTER  IMPLEMENTATION  OF  THE  NETWORK 

The  equivalent  electrical  network  representing  the  composite 
dielectric  may  be  conveniently  solved  using  computer  methods.  An 
important  prerequisite  to  obtaining  a  solution  is  the  development 
of  an  appropriate  scheme  for  identifying  both  pixel  grid  nodes 
and  the  internodal  capacitances.  In  the  scheme  adopted  here,  the 
pixel  grid  nodes  are  numbered  sequentially  starting  from  the 
ground  electrode  and  terminating  with  the  excitation  electrode. 
The  internodal  capacitances  are  identified  as  elements  in  a  two- 
dimensional  array.  The  numbering  scheme  is  illustrated  in 
Figure  6. 

The  node  numbering  scheme  is  not  arbitrary,  but  is  chosen  in 
such  a  fashion  as  to  maximize  the  stability  of  the  numerical 
solution.  By  beginning  at  the  ground  electrode  and  moving  along 
planes  of  increasing  potential,  we  can  minimize  truncation 
errors.  The  number  of  nodes  in  the  pixel  grid  depends  on  the 
lateral  boundary  conditions.  The  node  numbering  schemes  for  both 
'insulating'  and  'periodic'  boundary  conditions  are  illustrated 
in  Figure  7. 
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Figure  6.  Pixel  grid  nodes  (round  brackets)  and  internodal 
capacitance  (sqaure  brackets)  numbering  scheme  for  the  case  of 
'insulating'  boundary  conditions. 
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Figure  7.  Node  numbering  scheme  in  the  case  of  a)  'insulating' 
and  b)  'periodic'  boundary  conditions.  In  the  case  of 
'insulating'  boundary  conditions,  the  internodal  capacitances 
along  the  starred  (*)  paths  are  obtained  from  series  connection 
of  the  constituent  internodal  capacitances.  In  the  case  of 
'periodic'  boundary  conditions,  the  nodes  at  one  boundary 
effectively  'wrap  around'  to  the  other  boundary. 


INTERACTION  MATRIX 

Expansion  of  equations  8  and  9  over  the  entire  pixel  grid 
will  yield  a  set  of  simultaneous  equations  which  must  be  solved 
to  obtain  the  unknown  nodal  potentials.  Once  the  nodal 
potentials  are  known,  the  net  dielectric  permittivity  is 
calculated  as  the  quotient  of  the  excitation  displacement  current 
and  the  excitation  potential.  The  full  set  of  equations  takes 
the  form 


<  I  >  =  <  £  I  {  V  >  (10) 

where  (I)  is  a  column  vector  of  the  generalized  displacement 
currents,  (V)  is  a  column  vector  of  the  nodal  potentials,  and  (e) 
is  the  interaction  matrix  describing  the  composite  dielectric. 
The  interaction  matrix  is  in  general  a  sparse  symmetric  matrix, 
with  the  sparseness  arising  from  the  fact  that  only  nearest- 
neighbor  interaction  terms  are  non-zero.  Except  for  the  extrema, 
there  are  only  five  non-zero  terms  in  any  row,  hence  the  degree 
of  sparseness  increases  approximately  as  the  square  of  the  matrix 
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Figure  8.  Interaction  matrix  showing  the  location  of  non-zero 
terms  using  the  chosen  node  numbering  scheme. 


size.  The  location  of  the  non-zero  elements  in  the  matrix 
depends  upon  the  node  numbering  scheme  employed,  and  we  have 
chosen  a  node  numbering  scheme  such  that  the  interaction  matrix 
is  'banded'  with  the  non-zero  terms  clustered  about  the  main 
diagonal.  The  form  of  this  matrix  is  illustrated  in  Figure  8. 

Using  the  'banded'  interaction  matrix  as  shown  provides  a 
distinct  computational  speed  and  memory  size  advantage  over  the 
use  of  a  'non-banded'  matrix  with  arbitrary  non-zero  element 
locations.  The  execution  speed  in  the  'banded'  case  is 
proportional  to  the  square  of  the  matrix  size,  whereas  in  the 
'non-banded*  case  execution  speed  is  proportional  to  the  cube  of 
the  matrix  size.  Memory  size  is  conserved  using  the  'banded' 
form  since  only  a  small  segment  of  the  matrix  is  operated  upon  at 
any  time.  In  this  case,  the  size  of  the  interaction  matrix  which 
may  be  evaluated  is  limited  only  by  the  range  of  computer 
addresses  available  and  the  time  required  to  perform  the 
calculations. 
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VERIFICATION  OF  NUMERICAL  SOLUTION 


The  accuracy  of  the  numerical  solution  obtained  here  has 
been  tested  in  several  ways.  The  first  test  involved  calculation 
of  the  exponential  averaging  factor  a  for  the  known  cases  of 
stratified  layers  parallel  to  the  excitation  and  perpendicular  to 
the  excitation.  In  both  cases,  for  permittivities  e-,  and  e2  such 
that  z 2/ € 2  <  (limited  by  the  computer  implementation),  the 
results  of  the  numerical  solution  were  as  expected.  The  second 
test  involved  examining  the  stability  of  the  numerical  solution 
as  the  grid  size  was  changed  for  a  fixed  distribution  of 
constituent  materials.  No  variations  in  the  numerical  output 
were  detected  as  the  grid  size  was  varied  from  4x4  pixels  (26 
nodes)  to  40x40  pixels  (3161  nodes) . 

The  third  test  involved  comparing  numerical  solutions 
employing  different  implementations  of  node  numbering  and  matrix 
inversion  routines.  Essentially,  the  final  version  of  the 
numerical  analysis  code  grew  out  of  three  earlier,  less  efficient 
but  computationally  accurate  implementations.  The  first  code 
employed  a  node  numbering  scheme  which  started  at  the  center  of 
the  dielectric  and  spiraled  outward  to  the  boundaries.  This  code 
also  employed  the  matrix  inversion  routine  resident  in  the 
programming  language.  The  second  version  employed  the  same  node 
numbering  scheme  but  substituted  L-U  decomposition  for  matrix 
inversion.  The  third  version  also  used  the  same  node  numbering 
scheme,  put  employed  Gaussian  elimination  in  evaluating  the 
interaction  matrix.  All  three  versions  employed  back 
substitution  to  verify  the  numerical  solutions.  No  differences 
were  found  between  the  test  case  solutions  obtained  by  these 
versions  and  the  final  version. 

The  fourth  and  final  case  involved  comparing  the  current 
numerical  results  to  those  obtained  in  an  earlier  work  [3]. 
Again,  the  current  numerical  results  were  in  good  agreement  with 
the  earlier  work. 

PRELIMINARY  RESULTS 

The  computer  code  has  been  used  to  examine  the  variation  in 
exponential  averaging  factor  a  with  the  ratio  of  constituent 
permittivities  in  a  two-component  composite  dielectric.  Figures 
9  through  12  show  the  results  obtained  for  e,/e2  ratios  of  1.1, 
10,  100,  and  1000  respectively.  Each  figure  snows  the  value  of  a 
as  a  function  of  constituent  volume  ratio.  Each  data  point 
represents  the  calculated  behavior  of  a  capacitor  with  the 
indicated  volume  fractions  of  constituent  dielectrics.  For  each 
calculation,  the  spatial  distribution  of  the  constituent 
dielectrics  is  determined  by  a  random  number  generator,  the 
interaction  matrix  is  evaluated  to  determine  the  net  dielectric 
permittivity,  and  Equation  (1)  is  solved  iteratively  for  a.  One 
can  see  from  these  figures  that  the  behavior  of  the  exponential 
averaging  factor  depends  on  the  constituent  permittivity  ratio 
and  that  the  behavior  falls  into  three  general  categories: 
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Figure  9.  Exponential  averaging  factor  a  as  a  function  of 
constituent  volume  ratio  V-^/Vj  for  permittivity  ratio 
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Figure  10.  Exponential  averaging  factor  a  as  a  function  of 
constituent  volume  ratio  V-jyv2  for  permittivity  ratio  e2/el  =  10 • 
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Figure  11.  Exponential  averaging  factor  a  as  a  function  of 
constituent  volume  ratio  V-,/V2  for  permittivity  ratio 
€ 2/ e  =  100  . 
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Figure  12.  Exponential  averaging  factor  a  as  a  function  of 
constituent  volume  ratio  Vjyv2  for  permittivity  ratio 
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a)  f i/f 2  ~  1  -  the  exponential  averaging  factor  a  is 
essentially  constant  and  of  value  «0.25. 

b)  2  <  e1/e2  -  100  “  a  is  approximately  linearly  dependent 
on  the  constituent  volume  ratio. 

c)  e i/ e 2  >100  -  a  exhibits  a  nonlinear  constituent  volume 
ratio  dependence  arising  from  percolation  effects. 

Several  more  cases  will  need  to  be  examined  to  determine  more 
accurately  the  limits  of  these  behavioral  categories. 

FUTURE  WORK 

The  identification  of  the  exponential  averaging  factor  a  as 
falling  into  behavioral  categories  dependent  upon  the  constituent 
permittivity  ratio  is  of  considerable  engineering  significance. 
Once  the  behavioral  boundaries  are  established,  it  will  be 
possible  to  develop  simple  expressions  approximating  a  for  the 
various  behavioral  categories.  Future  work  in  this  area  should 
thus  include: 

a)  Evaluation  of  a  sufficient  number  of  test  cases  to 
determine  the  limits  of  the  behavioral  categories. 

b)  Development  of  engineering  tables  and  approximations  for 
practical  determination  of  a  in  engineering  applications. 

c)  Extension  to  the  case  of  complex  permittivity,  including 
artificial  dielectrics  such  as  dielectric-metal  composites. 

d)  Extension  to  3-dimensional  sample  geometries. 

e)  Analysis  of  the  nature  of  percolation  effects  observed 
for  large  constituent  permittivity  ratios,  in  particular, 
analysis  of  the  effects  of  constituent  grain  size. 

f)  Analysis  of  composites  with  three  or  more  components. 
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Appendix  I 

Two-Dimensional  Numeric  Analysis  Program 


10  !<<<<<<<  “0IEL_LR62O*  >>>>>>>> 

20  I*-#-*-#-#-*-*-*-#-*-*-*-#-*-*-#-*-* 

30  !  A  main  program  to  evaluate  a  2  dimensional  composite  dielectric 

40  1  response  for  a  pixel  network  of  capacitors. 

50  '  S.  R.  Wallin,  6/5/90 

60  I*-*-*-*-*-*-*-*-#-#-#-#-*-*-#-#-*-* 

70  PRINT  "  MEMORY  IS"  ;UAL(  SYSTEMS(  "MAILABLE  MEMORY"  )  )/8j  "<  reals  >" 

80  OPTION  BASE  1 

90  DATA  1,2,4,7,11,16,22,29,37  *  Dielectric  data  for  option 

100  COM  /Pass/  Relay  !  for  sharing  to  subs 

110  COM  /Pixel/  ChdrS£801 ,Ohdr*[80J .INTEGER  Lx tnt  ,Pix 1 ( 1 : 1 80 , 1 : 1 80 > 

120  DIM  Hp iv< 1 : 202  )  ,Hpr( 1 : 20302  )  !  dim  to  reasonable  size 

130  DIM  Hdr$C80]  !  available  string  for  headers 

140  DIM  Admt ( 0 : 3  )  ,Admsav( 0 : 3  )  1  neighbor  admittance  values 

150  INTEGER  Xt(0:3 )  ,Yt<0:3 )  !  neighbor  addresses 

160  COM  /Mem r/  Graf ( 1 : 256 , 1 i 4  )  ,Ahdr$C 803  ,Bhdr$I 80 1  ,  INTE6ER  Rep ,Kwd  (trials  mem 

170  i**#>  COM  areas  can  be  reaccessed  with  next  RUN  if  identical  name  &  sizes 

180  !»**>  nb.  ,  max  Lside  >-  .5  +  sqrt< .25  +  2*(max  dim  -  1)  > 

190  LET  Start-TIMEDATE 

200  INTEGER  Lside ,Kond  .Nodesz 

210  INTEGER  Ptrn.Nd.Ndl  ,Nd2  ,Xkin ,Ykin  ,Xcnt  ,Ycnt  ,Xaddr  .Yaddr  .Boxes  .Slant  .Sprss 

220  INTEGER  Qdrnt  ,Rptr  .Trans  , Pose ,6rpt  ,Tls  ,Sctr  ,Itmp  ,0ccp 

230  INTEGER  Hzmax  .Hxymax  .Hpremax .Hopped  .Hsteps  ,Hnd  ,Hcnt  ,Zcnt ,Znd 

240  INTEGER  Cnmb  .Qmat  ,Hnd1  ,Hnd2 ,Hlef t ,Hr ight 

250  DIM  Oiel( 1 :9),Frpx(0:9).Msd$I60] 

260  OIM  FI n$[ 60 ] 

270  (OIM  Hpiv< 1 : 128  )  ,Hpr< 1 : 8 1 92  )  !  set  to  max  physical  storage 

280  LET  Grpt»1  !  Initialize  the  data  storage  counter 

290  i, 

300  1  for  which  the  integer  variables  roles  are: 

310  1  Relay  -  an  available  common  pass  variable 

320  1  Lside  -  the  t  of  pixel  capacitor  elements  encounter  along  an  edge 

330  !  of  the  square  of  pixels 

340  I  Tls  ■  Lside  or  Lside/2  if  2x2  tiling 

350  1  Qdrnt  -  quadrant  pixel  array  expanding  switch,  0-off  &  1-on 

360  !  Px_tot  -  total  f  of  pixels  in  square  -  Lside*Lside 

370  •  Kond  -  the  boundary  condition  on  the  sides  of  the  overall  composite 

380  !  capacitor,  1>  insulating  sides  or  2)  periodic  or  sides  which 

390  1  wrap  around 

400  !  Nodesz  -  the  maximum  number  of  interaction  nodes  in  forming 

410  i  network,  with  0  as  the  ground  or  base  plate,  t  as  center 

420  1  node,  and  the  final  node  number  for  the  top  plate. 

430  I  Its  value  is:  L*L/2-L+2  +  (IF  Perlodic-1  )*( L/2- 1  ) . 

440  * 
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450  !  OVERALL  CAPACITOR  FROM  NOOES 

460  !  ! 

470  I  I 

480  !  mmmmmmmmmmmmmmmmmmmmm  top  n0£je  0  f  plate 

490  I  XXXXXXXXXX 

500!  XXXXXXXXX 

510  i  X  X  X  X  X  X  X  X  X  X 

520!  XXXXXXXXX 

530  !  XXXXXXXXXX  center  node  at  Midpoint 

540  i  XXXXXXXXX 

550  I  XXXXXXXXXX 

560!  XXXXXXXXX 

570!  XXXXXXXXXX 

580  I  mmmmmmmmmmmmmmmmmmmmm  base  node  OT  plate 

590  !  ' _ 

600  !  _ 

610  1  _  ...  X’s  represent  nodes 

620  ' 

630  I  Pixl<  )  ■  the  overlaying  matrix  representing  the  capacitor  pixels 

640  i  Osplc(  )  -  displacement  current  of  pixel  per  normalized  volt/meter 

650  !  Potnt<  )  ■  pixel  voltage  relative  to  one  volt  across  entire  sample 

660  !  Diel<  )  ■  dielectric  value  or  admittance  value  of  a  capacitor  pixel 

670  i  attached  to  addresses  represented  in  the  pixel  grid 

6€0  !  Frpx<  )  -  volume  fractions  associated  with  pixel  types 

690  !  Xt(  ),Yt(  )  -  neighbor  addresses 

700  !  Fln$  -  string  refering  to  a  filename.  Hdr®  -  80  chars 

7)0  •  AhdrS  ,Bhdr$  »  headers  of  80  chrs  for  Data  Title  &  ID  for  COM  /Mem r/ 

720  •  Chdr$,Dhdr$  «  headers  of  80  chrs  for  Pixel  Title  &  ID  for  COM  /Pixel/ 

730  !  Ptrn  -  choice  of  pixel  grid  filling  pattern 

740  !  Nd  »  a  single  number  label  for  a  node 

750  I  Ndl  ,Nd2  ■  refers  to  a  1st  node  &  a  2nd  node  0  specified  by  single 

760  !  node  numbers 

770  !  Xkin.Ykin  *  kinship  2D  address  of  a  node  number  i.e.  (x.y) 

780  *  Xcnt  .Ycnt  ■  step  counters  to  pixels  neighbouring  a  node  in  2D 

790  1  Xaddr.Yaddr  -  addresses  of  neighbouring  pixels  in  2D 

800  1  Boxes  •  total  concentric  boxes  fitting  within  pixel  grid  or 

810  !  number  of  2x2  cell  blocks  along  an  edge  on  pixel  grid 

820  1  Slant  ■  0  if  foward  slash  or  1  if  backslash  slanting  capacitor 

830  1  Rep  .Rptr  -  overall  number  of  repeats,  Kwd"t  of  data  storage  types 

840  i  6rpt  *  overall  plus  transpose  repeats  for  use  of  data  storage 

850  1  Sprss  »  Suppression  of  printout  details 

860  !  Trans  .Pose  ■  Pixel  transpose  selection 

870  *  Resp  "  Overall  dielectric  response  of  pixel  sample  along  E 

880  1  Tmp  .Tmpl  ,Tmp2  ,Vt 1  ,Vt2  "  reals  available  for  various  uses 

890  i  * 

300  PRINT 

910  PRINT  “  >  >  >  Happy  capacitor  composite  adventures  in  2  dimensions  <  <  <" 
920  PRINT  “  preformed  on  " i DATES!  TIMEDATE  ) ; 

930  PRINT  "  at  “ i T IMES< T IMEDATE  ) 

940  PRINT 

950 

960  l  The  hopper  reduction  subarray:  S.  Wallin,  July  1990 

970  !  .  .  >  large  symmetric  sparse  matrix  . 

980  •  !  1  ,1  ! 
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990  i 

!  9  1  _ 

1000  i 

!  2,1 

2  .2 

1010  i 

!  #2_ 

#3 

1020  i 

!  3.1 

3.2 

1030  i 

l_#4_ 

_#5_ 

The  hopper  subarray  moves 
down  as  pivoting  progresses 


\/  lg  sym  sparse  matrix 


1  The  "H“  prefix  is  mainly  used  to  denote  variable  use  in  hopper  program 
i  Hr...  -  row  address  of  array  element 
!  Hy..  -  column  address  of  array  element 
i  i.e.  (Hx..,Hy..)  specifies  a  coordinate  position 
*  Hz..  -  the  storage  number  associated  with  ( Hx . .  ,Hy . .  ) 
i  Hxymax  *  altitude  or  base  or  diag  f  elements  of  reduction  hopper 
i  Hzmax  -  the  total  #  of  elements  contained  within  reduction  hopper 
i  Hpremax  »  same  as  Hzmax  but  less  Hxymax  ( le  less  largest  row) 
i  Hsteps  -  extent  of  larger  interaction  matrix  upon  hopper  reduces 
I  Hleft.Hright  =  node  #s  of  insulated  sides  Pixel  grid 

1  Hpiv(  )  *  Pivoting  vector  of  node  reduction 
1  Hpr<  )  -  working  hopper  array  of  matrix  reduction 

Ittttttttttttttttft 

PRINT  ">Try  hopping  along  to  a  solution  of  sparse  matrices  at  "i 
DISP  “10  to  be:  0)default  1 )lab  3.5""  2>lab  hardisc  3.4)A.B  office 
INPUT  "5)user  defined”  ,Nd 
IF  Nd<0  THEN  STOP 


1250  IF  Nd-0  THEN  Msd*-'"' 

1260  IF  Nd-1  THEN  Msd*-“:CS80, 700.1" 

1270  IF  Nd-2  THEN  Msd*-" :CS80  ,700" 

1280  IF  Nd-3  THEN  Msd*-" :CS80 .703 ,0“ 

1290  IF  Nd-4  THEN  Msd*-" : CS80  ,703 . 1 “ 

1300  IF  Nd-5  THEN  INPUT  "Name  (completely)  storage? "  .Msd* 

1310  IF  MsdSO" "  THEN  PRINT  RPT*( "  "  ,50 )  ;  ”  storage ''  I  Msd* 

1320  PRINT  “  The  pattern  choices  are:“ 

1330  PRINT  “  0)  internal,  via  COM  /Mem r/‘ 

1340  PRINT  “  I  >  from  file  storage" 

1350  PRINT  "  2)  every  pixel  filled  by  user" 

1360  PRINT  “  3)  random  (i.e.  well  mixed)" 

1370  PRINT  “  4 >  by  slanted  fill  level” 

1380  PRINT  *  5)  with  an  circle  or  ellipse  of  which  can  be  tilted" 

1390  PRINT  "  6  )  with  strata" 

1400  PRINT  "  7)  concentric  boxes" 

1410  PRINT  “  8)  an  ellipse  with  host  &  inclusion  (2  components  only,": 

1420  PRINT  "  but  symm  wrt  1/2  vol  )" 

1430  INPUT  "Select  design  of  pixel  grid?  (see  above )"  ,Ptrn 
1440  IF  Ptrn<0  THEN  STOP 

1450  DISP  "  &  boundary  conditions?  1)  Insulative 

1460  INPUT  ”2)  Wrap  around  or  periodic "  ,Kond 

1470  IF  Kond<0  THEN  STOP 

1480  IF  Kond-0  THEN  Kond-1 

1490  IF  Kond>2  THEN  Kond- 1 +BI T( Kond+ 1  .0 > 

1500  INPUT  “Use  2x2  tiles  on  pixel  grid?  0 )  No  1)  Yes”,Tls 
1510  IF  Tls<0  THEN  STOP 
1520  LET  T1 s-1 +( Tls-1  ) 


17 


1530  INPUT  "Quad  fold  symmetry  expansion  of  pixel  grid7  0 >  No  1  )  Yes"  .Qdrnt 
1540  IF  Qdrnt <0  THEN  STOP 
1550  LET  QdrntH  +  <  Qdrnt"!  ) 

1560  SELECT  Ptrn 

1570  CASE  -0 

1580  IF  Lxtnt <2  THEN 

1590  PRINT  “  Are  Pixels  there  in  memory7  ..idled  ..start  again" 

1600  STOP 

1610  END  IF 

1620  IF  Lxtnt/181  THEN  PRINT  “  ..  there  may  be  too  many  Pixels" 

1630  REDIH  Pixl ( 1 :Lxtnt  ,  1 :Lxtnt  ) 

1640  LET  Lside-Lxtnt»Qdrnt»Tls 

1650  PRINT  "  From  internal  memory  via  COM,  Pixels" iLxtnt ; “x“ ;Lxtnt 5 “  ,  title, 
1660  IF  Chdr$<> " "  THEN  PRINT  Chdr# 

1670  IF  DhdrSO "  “  THEN  PRINT  OhdrS 

1680  CASE  ”1  1  6et  pixels  from  file 

1690  INPUT  "  Enclose  (in  . s)  file  name  to  contain  pixel  pattern7" , Fin# 

1700  IF  Fin*-""  THEN  STOP 

1710  IF  POStFlnS  :  “  )-0  THEN  Fln$-Fln#8,Msd$ 

1720  DISP  "  File  named  . iFln*i“*“  < C “ ;LEN< Fin# > : "  ]  characters  ) “  ; 

1730  OISP  *  is  being  read  from  storage" 

1740  ASSIGN  @P i x src  TO  Fin# ; FORMAT  OFF 

1750  ENTER  @Pixsrc iChdr#  .DhdrSiLxtnt  !  NB  header  assigned  length  of  80 

1760  PRINT  "  Pixels  contained  in  file  ""“;Fln$; . ,  entitled  with” 

1770  PRINT  Chdr# 

1780  PRINT  Dhdr# 

1790  REDIM  Pix 1 ( 1 :Lxtnt  ,  1 iLxtnt  )  I  read  initial  Pixl(*>  array 

1800  ENTER  @Pixsrc  ?  Pi  x  K  *  )  1  retrace  pixels  from  file 

1810  ASSIGN  QPixsrc  TO  *  I  close  file 

1820  LET  Lside-Lxtnt*Qdrnt*Tls  i  actual  Pixel  side  anticipated 

1830  PRINT 

1840  CASE  ELSE  1  Generate  pixels 

1850  DISP  "How  big  a  capacitor  pixel  grid  in  elements/side7  “; 

1860  INPUT  "(even  t,  max  '180  int  addr  lmt >" ,Lside 
1870  IF  Lside<0  THEN  STOP 

1880  IF  Ls i de<>SHIFT( SHIFT< Lside  ,  1  > ,  —  1  )  THEN 

1890  PRINT  “  Odd" ;Lside; "Pixel  length  changed  to  even"; 

1900  LET  Lside-SHIFT(SHIFT(Lside,1  >  .- 1  ) 

1910  PRINT  Lside 

1920  END  IF 

1930  IF  Lside-0  THEN  Lside*2 

1940  LET  Lxtnt"Lside  1  initial  Pixel  side  length 

1950  LET  Lside“Lside«Qdrnt«Tls  1  Pixel  side  length  anticipated 

1960  IF  Lside>182  THEN  PRINT  "  ...  near  integer  addressing  limit" 

1970  END  SELECT  i  end  Ptrn  test 

1980  PRINT 

'990  IF  Tls«2  OR  Qdrnt»2  THEN  PRINT  "  Pixels  now  measures" iLside; “x” sLside 
2000  PRINT  “  Pat  tern* "; Ptrn; ")“; Ls ide x " ;Ls l de ; 

2010  IF  Tls-1  THEN  PRINT  "pixels,"; 

2020  IF  Tls-2  THEN  PRINT  "tiled  pixels,"; 

2030  IF  Kond-1  THEN  PRINT  "  insulated  or  H“D“"  field  parallel  to  edge." 

2040  IF  Kond*2  THEN  PRINT  "  periodic  or  voltage  wrapping  around  at  edges." 

2050  LET  Nodesz-SHIFT(Lside*Lside  ,1  )-Lside  +  2  +  ( Kond-2  )•( SHIFT! Lside  ,  1  )-1  > 

2060  lALLOCATE  REAL  Osp lc( Ls ide  ,Ls l de  )  ,Po t nt < Ls i de  ,Ls l de  )  •  to  be  programmed 
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2070  !*»*>  Initializing 

2080  MAT  Oiel-  (0) 

2090  DISP  "Dielectric  sources?  “s 

2100  OISP  “1  >  data  l nput ( programmed  )  2)  keyboard  "i 

2110  INPUT  ”3)  by  progress  ion"  ,Nd 

2120  IF  Nd-0  THEN  Nd-2 

2130  SELECT  Nd 

2140  CASE  -1 

2150  FOR  Ndl-1  TO  9 

2160  READ  Diel( Ndl  > 

2170  NEXT  Ndl 

2180  CASE  -2 

2190  FOR  Ndl-1  TO  9 

2200  DISP  "Give  dielectric  value  a t “ j Nd 1 ; 

2210  INPUT  “?  (or  enter  negative  if  to  cease  )"  ,Diel( Ndl  ) 

2220  IF  DieKNdl  K0  THEN 

2230  LET  DieKNdl  )-0 

2240  LET  Nd 1 -9 

2250  ELSE 

2260  PRINT  "diell"; Ndl »"]-"» PR0UND( Diel ( Ndl  >  ,-4  > ;  “  ,  “ i 

2270  END  IF 

2230  NEXT  Ndl 

2290  CASE  -3 

2300  INPUT  "Dielectric  value  of  pixel  type  " "[ 1 ] " ,Tmp 

2310  OISP  “Multiplier  of  progression  for  each  succeeding  value  “  j 

2320  INPUT  "to  fill  C  2  1  ,  C  3 1 _ C91  7“  .Trnpl 

2330  FOR  Ndl-1  TO  9 

2340  DieKNdl  >-Tmp 

2350  Tmp-Tmp*Tmp  1 

2360  PRINT  “die It “ (Ndl i  “  l-“ ;PR0UND( DieK Ndl  )  ,-4  K  "  ,  “  5 

2370  NEXT  Ndl 

2380  CASE  ELSE 

2390  STOP 

2400  END  SELECT 

2410  PRINT 

2420  LET  Sprss-1 

2430  IF  Nodesz<32  THEN 

2440  INPUT  “Surpress  screen  listing  details,  0 )  No  1)  Yes7"  ,Sprss 

2450  IF  Sprss<0  THEN  STOP 

2460  END  IF 

2470  INPUT  “Any  overall  repeat s7"  , Rep 
2480  IF  Rep<0  THEN  STOP 
2490  IF  Rep-0  THEN  LET  Rep-1 

2500  INPUT  "Oesire  transpose  of  pixel  grid7  0 )  No  1)  Yes", Pose 
2510  IF  Pose<0  THEN  STOP 
2520  LET  Pose- 1 +BI T<  Pose  ,0  ) 

2530  REM  "  Solution  acheived  by  a  sparse  matrix  reduced  pivoting  technique" 
2540  !*#*>>  Overall  repetition,  may  require  additional  editing 

2550  LET  Kuid-4  1  user  has  selected  to  program  for  4  data  columns 

2560  REOIM  Graft  1 :Rep  ,1 :Kwd) 

2570  FOR  Rptr-1  TO  Rep 

2580  LET  Relay-Rptr  '7/27/90 

2590  MAT  Frpx-  <0> 

2600  !**>  if  Ptrn-0  internal  or  Ptrn-1  then  Pixels  read  from  file 
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2610  IF  Ptrn-2  THEN  CALL  Pixl2d_fill 

2620  IF  Ptrn-3  THEN  CALL  Pixl2d_rand 

2630  IF  Ptrn-4  THEN  CALL  Pixl2d_tilt 

2640  IF  Ptrn-5  THEN  CALL  Pixl2d_ellP3 

2650  IF  Ptrn-6  THEN  CALL  Pixl2d_strat 

2660  IF  Ptrn-7  THEN  CALL  PixI2d_cbox 

2670  IF  Ptrn-8  THEN  CALL  Pixl2d_ellp2 

2680  IF  Tls-2  THEN 

2690  LET  Xkin**SIZE(Pixl  ,  1  >  *  redimensioning 

2700  LET  Ykin-SIZE(Pixi  ,2  ) 

2710  LET  Ndl-Xkin+Ykin 

2720  LET  Nd2-SHIFT(Xkin*Ykin,-2>  '  4«Xkin*Ykin 

2730  REDIM  Pixl( 1:1,1 :Nd2  ) 

2740  FOR  Xcnt*<  Xk in- 1  )  TO  0  STEP  -l 

2750  FOR  Ycnt-Ykin  TO  1  STEP  -1 

2760  LET  Pixl( 1  ,Xcnt*2*Yk m+Ycnt  )*Pixl( 1  ,Xcnt*Ykin+Ycnt  ) 

2770  NEXT  Ycnt 

2780  NEXT  Xcnt 

2790  REDIM  Pixl< 1 :Ndl  ,1 : N d 1  )  1  set  new  array  diner 

2800  FOR  Xcnt-Xkin  TO  1  STEP  -1  '  tiling  2x2 

2810  FOR  Ycnt-Ykin  TO  1  STEP  -1 

2820  LET  Itmp*Pix K Xcnt  ,Ycnt  ) 

2830  LET  Xaddr=*SHIFT(  Xcnt  1  )  *  effective  2*  op 

2840  LET  Yaddr-SHIFK Ycnt  ,-1  ) 

2850  LET  PixKXaddr  .Yaddr  >-Itmp 

2860  LET  PixKXaddr — I  , Yaddr  )-Itmp 

2870  LET  PixKXaddr  .Yaddr-1  )-Itnp 

2880  LET  Pixl(Xaddr-l  , Yaddr-1  )-Itmp 

2890  NEXT  Ycnt 

2900  NEXT  Xcnt 

2910  END  IF 

2920  IF  Qdrnt-2  THEN 

2930  LET  Nd2-SIZE(Pixl  ,1  )+SIZE(Pixl  ,2) 

2940  LET  Ndl “SHIFT ( Nd2  .1  ) 

2950  REDIM  P i x K 1 : Nd2  , 1 : Nd2  )  f  redim  to  dble  quad  duplic 

2960  FOR  Xcnt "Nd 1  TO  1  STEP  -I 

2970  LET  Xaddr*Nd2+ 1 -Xcnt  1  quad  complement  X  counter 

2980  LET  Xkin-SHIFT< Xcnt+1  ,1  >  '  effectively  DIO  2  op 

2990  LET  Yk in*BIT( Xcnt  +  1  ,0  )  •  effectively  odd<*>even  op 

3000  FOR  Ycnt* 1  TO  Ndl 

3010  LET  Yaddr*Nd2+ 1 -Ycnt  1  quad  complement  Y  counter 

3020  LET  Itmp*Pixl(Xkin  ,  Ycnt +Nd 1  * Yk in  ) 

3030  LET  Pix 1( Xcnt  , Ycnt  )-Itmp  '  Itmp  takes  care  of  redim  elements 

3040  LET  Pi x K Xaddr  , Ycnt  )*I tmp 

3050  LET  P i x K Xcnt  , Yaddr  )- 1 tmp 

3060  LET  Pix K Xaddr  , Yaddr  )*I tmp 

3070  NEXT  Ycnt 

3080  NEXT  Xcnt 

3090  END  IF 

3100  LET  Ls ide-SI ZE< Pix  1  , 1  )  t  update  Pixl  extent  along  edge 

3110  LET  Px_tot“Lside»Lside 

3120  LET  Boxes-SHIFT(  Ls i de , 1  ) 

3130  IF  Boxes<1  THEN  PRINT  "WARNING'  may  be  too  small  of  a  pixel  grid" 

3140  LET  Nodesz*SHIFT<Lside*Lside,l  >-Ls i de  +  2  +  < Kond-2 )* ( SHIFK Ls i de . 1  >-!  ) 
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3150 

3160 

3170 

3180 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3390 

3400 

3410 

3420 

3430 

3440 

3450 

3460 

3470 

3460 

3490 

3500 

3510 

3520 

3530 

3540 

3550 

3560 

3570 

3580 

3590 

3600 

3610 

3620 

3630 

3640 

3650 

3660 

3670 

3680 


i***>  lranpose  of  Pixel  grid 
LET  Trans-Pose  *  if  loop  then  use  next  line 
'FOR  Trans-1  TO  Pose 

IF  Trans-2  THEN  1  Tranpose 

PRINT  "  TRANSPOSING" 

FOR  Xcnt-1  TO  Lside 

FOR  Ycnt*(Xcnt+1 )  TO  Lside 

LET  Itmp-Pixl< Ycnt ,Xcnt >  1  swap  Xcoordinate<->Ycoordinate 

LET  PixKYcnt  ,Xcnt  )-Pixl(Xcnt  ,Ycnt  ) 

LET  Pix  1  (  Xcnt  ,Ycnt  )-I  trip 
NEXT  Ycnt 
NEXT  Xcnt 
ENO  IF 

1**0  Evaluation  of  pixel  type  volume  fractions 

IF  NOT  (Sprss)  OR  Lside<81  THEN  PRINT  “  P l xe  1  s “ i Ls i de ; " x " ; Ls i de 
FOR  Xcnt=*Lside  TO  1  STEP  -1 
FOR  Ycnt-1  TO  Lside 


Frpx(Pixl(Xcnt , Yen t  ))-Frpx(Pixl(Xcnt , Yen t ) )  +  1 

IF  NOT  (Sprss)  OR  Lside<40  THEN  PRINT  “  VAL$( P l x 1 ( Xcnt , Ycnt  )> i 
NEXT  Ycnt 

IF  NOT  (Sprss)  OR  Lside<40  THEN  PRINT 
NEXT  Xcnt 

HAT  Frpx-  Frpx / ( Px_t ot  ) 

PRINT  “  Volume  Xs :  “ i 
FOR  Nd-1  TO  9 

IF  Frpx (  Nd  ) <  >0  THEN  PRINT  PROUND(  I  00*Frpx(  Nd  )  ,  -  1  );"%->["; Nd ; 

NEXT  Nd 
PRINT 

IF  Frpx ( 0  >< >0  THEN  PRINT  "WARNING'  check  pixels" 

OISP  “  ..  wai t " i Rp t r ; "of " ( Rep i ” . .  solving  node  interact  matrix,'; 

OISP  Nodes. ; "by" ;Nodesz ;“ from  time  "  ;  TIME$( TIMEDATE > 

LET  Tmp-TINEDATE  1  Benchmarker 

'*•*>  HOP  technique  of  sparse  matrix  reduction 

PRINT  "  Solving  INTERACTION  matr i x  ”; Nodes; x Nodes v i a  hopper” 

LET  Hvymax-Lside  '  size  of  Hopper  'or  matrix  reduction 

LEt  Hzna < - (  1 . 0+Hx ymax  )  *Hx ymax  DIO  2  '  memory  reduction  capabilities 

lET  Hsteps-Nodesz 

lE T  Hpremax -Hzmax-Hx ymax 

REDIH  Hpiv(  1  :Hxymax  )  ,Hpr( I :Hzmax ) 

lET  Hieft-1  '  markers  of  node  8s  on  left  insi  BC 


LET  Hright-Boxes  1  markers  of  node  8s  on  right  insi  BC 

MAT  Hpiv-  (0) 

MAT  Hpr-  ( 0  ) 

LET  Znd-0  1  to  previous  row  tn-diag  accum 

FOR  Hnd-1  TO  Hxymax  i  Filling  hopper  work  array 

IF  NOT  (Sprss)  THEN  PRINT  “  node's  ‘ ; VAL$( Hnd > i "  neighbors  are"; 

FOR  3ctr»0  TO  3  '  Diagonal  or  self  interai  t  terms 

IF  Kond-1  THEN  CALL  C vnd i ( Hnd  ,Ls i de  ,Sc t r  ,X t ( Sc t r  )  ,  Y t ( Sc t r  )  ) 

IF  Kond-2  THEN  CALL  C vndp < Hnd ,L s i de , Sc t r  , X t ( Sc t r  )  ,  Y 1 1 Sc t r  ) ) 

LET  AdmtiSctr  /-Diel(Pixl(Xt(Sctr  )  , Y t ( Sc tr  )  )  ) 

IF  Kond-1  THEN  >  adj  for  insi  BC  on  Pixel  grid 

IF  Hnd-Hieft  THEN  '  test  if  left  side  node 

IF  NOT  <  BIT<  5ctr  ,0  )  )  THEN 

IF  BIT(Sctr,1 )  THEN  •  make  series  admit  combo 
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3690 
3700 
3710 
3720 
3730 
3740 
3750 
3760 
3770 
3780 
3790 
3800 
3810 
3820 
3830 
3840 
3850 
3860 
3870 
3880 
3890 
3900 
3910 
3920 
3930 
3940 
3950 
3960 
3970 
3980 
2990 
4000 
4010 
4020 
4030 
4040 
4050 
4060 
4070 
4080 
4090 
4100 
41  10 
4120 
4130 
4140 
4  150 
4160 
4170 
4180 
4  190 
4200 
4210 
4220 


IF  X t  (  2  XLside  THEN 
LET  Admsav(0)*Diel(pixl(Xt(2  )+1  ,1  )) 

LET  Admt ( 2 )*Admt  <  2 >*Admsav( 0  )/ ( Admt ( 2  )+Admsav<  0  )  ) 

LET  Admsav( 2  )“Admt < 2 >  !  save  for  later  use 
END  IF 

ELSE  !  reuse  former  upper  admit  combo 

IF  Hnd>1  THEN  LET  Admt ( 0 >°Admsav( 2 >  1  pass  by  on  1st  node 
END  IF  '  end  of  Sctr-0,2  test 

END  IF  1  end  of  msl  BC  Pixel  test 

END  IF  1  end  of  Hnd*Hleft  test 

IF  Hnd-Hright  THEN  1  test  if  right  side  node 

IF  B I T < Sctr  ,0  )  THEN 

IF  BIT ( Sc tr ,  1  >* 1  THEN  1  make  series  admit  combo 


IF  X t (  3  XLside  THEN 

LET  Admsav< 1  )-Diel < Pixl ( Xt < 3 >+ 1  ,Ls i de >  )  1  cmbn  above  ngbr 
LET  Admt(3)"Admt(3)*Admsav(1>/<Admt(3>+Admsav(1  1  ) 

LET  Admsa v( 3  )”Admt ( 3 )  <  save  for  later  use 


END  IF 
ELSE 

IF  Hnd>Boxes 
END  IF 
END  IF 
END  IF 
END  IF 


1  reuse  former  admit  combo 
THEN  LET  Admt ( I  )-Admsav( 3  )  1  Skip  1st  pass 
1  end  of  Sctr-0  test 
i  end  of  msl  BC  Pixel  test 
!  end  of  Hnd*Hnght  test 
i  end  if  for  Kond»1  test 


LET  Hpr ( Znd+Hnd  >”Hpr( Znd+Hnd )+Admt  <  Sctr  ) 

IF  NOT  ( Sprss  )  THEN  PRINT  "  < “ ! UAL$< X t ( Sc tr  )  )  i “  , " ; UAL$< V t ( Sc t r >  ) ; " 
NEXT  Sctr 
SELECT  Kond 
CASE  - i 

LET  Hnd1-FNNi(Xt<0)  ,Yt<0)  .Lside.l  > 

LET  Hnd2-FNNi<  Xt( 1  )  ,Yt( 1  )  , Lside  ,1  > 

CASE  -2 

LET  Hndl -FNNp< Xt ( 0 >  ,Yt<0  ) , Lside  ,  1  > 

LET  Hnd2=*FNNp<  X  t  (  1  )  ,Yt(  1  )  .Lside  ,1  > 

CASE  ELSE 

PRINT  "  out  of  bounds,  boundary  condition,  in  HOPper" 

END  SELECT 

IF  NOT  (Sprss)  THEN  PRINT  “  ui/  lower  nodes  “  i  Hnd  1  t Hnd2 
IF  Hnd1<0  OR  Hnd2<0  THEN  PRINT  "  Warning  node  addresses7” 

IF  Hndl < Hnd  ANO  Hnd1>0  THEN 

LET  Hpr(Znd+Hnd!  )“Hpr( Znd+Hnd 1  >-Admt<0) 

END  IF 

IF  Hnd2<Hnd  AND  Hnd2>0  THEN 


LET  Hpr( Znd  +  Hnd2  )”Hpr( Znd  +  Hnd2 )-Admt(  1  ) 

END  IF 

IF  Kond"1  THEN  •  increment  msl  BC  left  &  right  node 

IF  Hnd-Hleft  THEN  LET  Hlef t-Hlef t+Lside- 1 


IF  Hnd’Hr  lght  THEN  LET  Hr l ght »Hr l ght +L s l de- 1 
ENO  IF 

LET  Znd-Znd+Hnd  1  loop  count  accumulator 

NEXT  Hnd  1  end  of  set  up  of  work  matrix 

LET  Zcnt“0 


IF  NOT  (Sprss)  THEN 
FOR  Xcnt“1  TO  Hxymax 


i  printout  HOPper 


)“ 


f 
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4230 

FOR  Ycnt-1  TO  Xcnt 

4240 

PRINT  PROUND<  Hpr( Zcnt+Ycnt ) 

.-3). 

4250 

NEXT  Ycnt 

4260 

PRINT 

4270 

LET  Zcnt*Zcnt  •>  Xcnt 

4280 

NEXT  Xcnt 

4280 

END  IF 

4200 

FOR  Hopped*)  TO  Hatepa-) 

1  Let  the  pivoting  begin,  &  drip  dry 

43)0 

LET  Hnd*Hxymax+Hopped 

1  count  of  oncoming  node  number 

4320 

LET  Hp i v (  t  )* 1 

1  normalize  to  1  at  elment  pivot  vectr 

4330 

LET  Hnrm*l/Hpr( 1 ) 

1  normalizing  multiplier 

4340 

LET  Zcnt-2 

1  convert  array  atorage  for  )at  colmn 

4350 

FOR  Hcnt*2  TO  Hrynax 

1  set  pivot  vector 

4360 

LET  Hp  i  v< Hcnt  )*Hpr( Zcnt  )*Hnrm 

4370 

LET  Zcnt*Zcnt+Hcnt 

4380 

NEXT  Hcnt 

4390 

'•**>  one  can  output  pivot  here 

for  backaub  later 

4400 

IF  NOT  (Spras)  THEN 

4410 

PRINT  "  At  reduct  ion" iHoppedi 

"the  p i vot a  are : “ 

4420 

FOR  Hcnt*Hxymax  TO  1  STEP  -1 

4430 

PRINT  PROUND( Hpiv( Hcnt  )  ,-4 >i 

4440 

NEXT  Hcnt 

4450 

PRINT 

4460 

END  IF 

4470 

LET  Zcnt*) 

initialize  filling  counter 

4480 

FOR  Xcnt-2  TO  Hxymax 

heart  of  pivoting 

4490 

IF  Hp i  v(  Xcnt  K>0  THEN 

aparaeneaa  efficiency  *0,  no-op 

4500 

FOR  Ycnt*2  TO  Xcnt 

adj  each  array  row  with  pivot  vectr 

4510 

IF  Hpiv< Ycnt  )< >0  THEN 

aparaeneaa  efficiency  *0,  no-op 

4520 

LET  Hpr( Zcnt  +  Ycnt  >-Hpr<  Zcnt  +  Ycnt  )-Hpr( Zcnt+1  )*Hpiv<  Ycnt  > 

4530 

END  IF 

4540 

NEXT  Ycnt 

4550 

END  IF 

4560 

LET  Zcnt*Zcnt+Xcnt 

4570 

NEXT  Xcnt 

4580 

LET  Zcnt*0 

initialize  filling  counter  lower 

4590 

FOR  Xcnt*l  TO  Hxymax-1 

4600 

FOR  Ycnt-I  TO  Xcnt 

hopping  along  for  ahake  up 

4610 

LET  Hpr< Zcnt  +  Ycnt  )*Hpr( Zcnt  +  Ycnt  +  )+  Xcnt  > 

4620 

NEXT  Ycnt 

4630 

LET  Zcnt*Zcnt+Xcnt 

4640 

NEXT  Xcnt 

4650 

FOR  Ycnt*1  TO  H> y»a» 

feed  hopper,  clear  laat  row 

4660 

LET  Hpr < Hpremax+Ycnt )*0 

4670 

NEXT  Ycnt 

4680 

SELECT  Hnd 

4690 

CASE  YHatepa 

feed  unleaa  over  lg  array  extent 

4700 

IF  NOT  (Spraa)  THEN  PRINT  "  node'a  " i VAL$< Hnd  ) ; “  neighbora  are"  ! 

47)0 

FOR  Sctr-0  TO  3 

with  diagonal  or  aelf  interact  terma 

4720 

IF  Kond*)  THEN  CALL  Cvndi < Hnd  ,Ls l de  ,Sc t r  ,X t ( Sc tr  )  , Yt ( Sctr  )  ) 

4730 

IF  Kond*2  THEN  CALL  Cvndp< Hnd  ,La l de  ,Sc tr  ,X t ( Sc tr >  , Y t ( Sc tr )  ) 

4740 

LET  Adm  t  ( Sctr  )*Diel<  Pixl( Xt<  Sctr  >  ,  Y  t  <  Sctr  )  )  ) 

4750 

IF  Kond*)  THEN  t 

adj  for  mal  BC  on  Pixel  grid 

4760 

IF  Hnd-Hleft  THEN  i 

teat  if  left  aide  node 
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IF  NOT  ( BIT( Sctr  ,0  ) )  THEN 

IF  BIT(Sctr.l)  THEN  !  make  series  admit  combo 
IF  X  t  <  2  XLsioe  THEN 


LET  Admsav( 0  )“Diel ( Pixl <  Xt<  2  )+1  ,1  >  ) 

LET  Admt(2)“Admt(2>*Admsav(0)/(Admt(2)+Admsav(0)> 

LET  Admsav( 2  )“Admt ( 2  )  1  save  for  later  use 
END  IF 

ELSE  !  reuse  former  upper  admit  combo 

IF  Hnd>1  THEN  LET  Admt ( 0 >“Admsav< 2  ) 1  pass  by  on  1st  node 
END  IF  !  end  of  Sctr“0,2  test 

END  IF  !  end  of  insl  BC  Pixel  test 

END  IF  !  end  of  Hnd"Hleft  test 

IF  Hnd“Hright  THEN  <  test  if  riQht  side  node 

IF  BIT<  Sctr  ,0  )  THEN 

IF  BIT(Sctr,l  )*1  THEN  !  make  series  admit  combo 
IF  Xt( 3  XLside  THEN 

LET  Admsav( 1  )”0iel < Pix 1 ( Xt ( 3  )+1  ,Lside  )  ) 1  cmbn  above  ngbr 
LET  Admt ( 3 )“Admt ( 3  )*Admsav ( 1  )/(Admt(3)  +  Admsav(  1  )) 

LET  Admsavl 3 >“Admt ( 3  ) 1  save  for  later  use 
END  IF 

ELSE  1  reuse  former  admit  combo 


IF  Hnd>8oxes  THEN  LET  Admt(1  )“Admsav<3>1  Skip 


END  IF 
END  IF 
END  IF 
END  IF 


*  end  of  Sctr=*0  test 
!  end  of  insl  BC  Pixel  test 
1  end  of  Hnd-Hright  test 
1  end  if  for  Kond=*1  test 


1st  pass 


LET  Hpr <  Hzmax  )-Hpr( Hzmax  )+Admt ( Sctr  ) 

IF  NOT  ( Sprss  )  THEN  PRINT  "  < “ t  OALS  ( Xt  <  Sc tr  ) ) i " , " t  VALS( Yt ( Sc tr  >  ) ; 
NEXT  Sctr 

SELECT  Kond  !  with  off  diagonal  terms 

CASE  -1 

LET  Hnd1-FNNi(Xt(0),Yt(0),Lside,1  > 

LET  Hnd2*FNNi ( Xt( 1  )  ,Yt( 1  ) .Lside .1  > 

CASE  -2 


LET  Hnd I -FNNp( Xt<0),Yt<0),Lside,l  ) 

LET  Hnd2-FNNp( Xt( 1  )  , Y t ( 1  ), Lside  ,1  ) 

END  SELECT  '  to  SELECT  Kond 


IF  NOT  (Sprss)  THEN  PRINT  “  w/  lower  nodes"  jHndl  ;Hnd2 
IF  Hnd1<0  OR  Hnd2<0  THEN  PRINT  ”  Warning  node  <“>  addresses7" 
IF  HndKHnd  AND  Hnd1>0  THEN 
LET  Hndl-Hndl +Hzmax-Hnd 
LET  HpHHndl  )-Hpr(Hnd1  )-Admt(0> 

END  IF 

IF  Hnd2<Hnd  AND  Hnd2>0  THEN 
LET  Hnd2*Hnd2+Hzmax-Hnd 
LET  Hpr< Hnd2  )“Hpr( Hnd2 >_Admt ( 1  > 

END  IF 

CASE  “Hsteps  1  top  node  is  contact  node 

FOR  Ycnt»1  TO  Lside  1  diagonal 

Hpr( Hzmax )*Hpr ( Hzmax  )+Oiel(Pixl(Lside  ,Ycnt  >  ) 

LET  Hndl *Hzmax-Boxes+SHIFT( Ycnt- 1  ,1  )  1  off  diagonal  node  f 

LET  Hpr<  Hnd I  '“Hpr ( Hndl  )-Diel(Pixl(Lside,Ycnt)) 

NEXT  Ycnt 

END  SELECT  I  to  SELECT  Hnd 
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IF  NOT  (Sprss  OR  Hnd>Hsteps >  THEN 

FOR  Ycnt*1  TO  Hxymax  l  cont  printout  of  HOPper 

PRINT  PRQUNO( Hpr<  Hpremax+Ycnt  )  ,-3 ) , 

NEXT  Ycnt 
PRINT 
END  IF 

IF  Kond»1  THEN  1  increment  msl  BC  left  &  right  node  # 

IF  Hnd*Hleft  THEN  LET  Hleft*Hlef t+Lside-1 
IF  Hnd*Hnght  THEN  LET  Hr  lght  “Hr  lght  +Ls  l de- 1 
END  IF 
NEXT  Hopped 

IF  NOT  (Sprss)  THEN  PRINT  "  Hopper  funnels  down  to“;Hpr(1  > 

>***>  Hpr( 1  )  contains  the  end  of  the  interaction  reduction 
•  LET  Hp i v < Hs t ep s  )* 1 /Hpr ( 1  )  1  backsubst ltute  for  solution  vector 

1  FOR  Hcnt*( Hsteps-1 )  TO  1  STEP 
i  LET  Hpiv( Hcnt )*0 

1  FOR  Zcnt*Hcnt  TO  Hsteps 

!  LET  Hpiv(Hcnt  )*Hp 1 v( Hcnt  )-Pivotstorg(Hcnt  ,1+Zcnt-Hcnt )*Hpiv(Zcnt  > 

i  NEXT  Zcnt 

!  NEXT  Hcnt 

PRINT  "  ...  at  "  TIME$(  TIMEDATE  )  i "  inversion  excution  “  j 

PRINT  "time  t ook PR0UND< T IMEDATE-Tmp 1  ) 5 " seconds " 

DISP 

LET  Resp-Hpr( 1 )  !  principal  diel  resp 

'*•*>  Pixl  displacement  field/current  &  potentials 
1  FOR  Ycnt-1  TO  Lside 
1  FOR  Xcnt*l  TO  Lside 

1  Slant*! Xcnt+Ycnt )  MOO  2 

•  Xaddr-Xcnt-Boxes-Slant  1  (x,y)  of  node  upper  to  pixel 

1  Yaddr*Ycnt -Boxes 

1  CALL  Xy_to_node( Ndl  .Xaddr  ,Yaddr  ,Lside  ,Kond  ) 

1  Xaddr»Xcnt-Boxes+Slant-l  1  (x,y)  of  node  lower  to  pixel 

1  Yaddr*Ycnt-Boxes~ 1 

1  CALL  Xy_to_node( Nd2 .Xaddr ,Yaddr ,Ls ide ,Kond ) 

1  IF  Nd I < >Nd2  ANO  Nd1>0  AND  Nd2>0  THEN 

1  Dsplc(Xcnt,Ycnt)»0iel(Pixl(Xcnt,Ycnt))*(Hpiv(Nd1  )-Hpiv(Nd2)) 

1  Potnt ( Xcnt  , Ycnt >-<Hpiv<Nd1  >+Hpiv(Nd2  )  >/2 

1  END  IF 

'  IF  Nd2-0  AND  Nd1>0  THEN 

1  OsplclXcnt  ,Ycnt  )*Diel(Pixl(Xcnt  ,Ycnt  ))*Hpiv(Nd1  ) 

1  Potnt! Xcnt  .Ycnt  >-Hpiv(Nd1  >/2 

1  END  IF 

'  NEXT  Xcnt 

'  NEXT  Ycnt 

!*•*>  additional  modification  of  Potential  &  Displacement  array  fields 
1  IF  Kond-I  THEN 

1  FOR  Ycnt-2  TO  (Lside-2)  STEP  2 

1  FOR  Nd— 1  TO  I  STEP  2 

1  LET  Xcnt*<  Nd+!  )*Boxes  +  <  Nd— I  >  1  (Xcnt, Ycnt)  refer  to  pixel 

1  LET  Xaddr-Nd*( Boxes-1  )  1  ( Xaddr  ,Yaddr  )  refer  to  node 

1  LET  Yaddr-Ycnt-Boxes 

1  CALL  Xy_to_node( Ndl  , Xaddr  ,Yaddr+1  ,Ls 1 de  ,Kond  ) !  node  f  upper 

1  CALL  Xy_to_node( Nd2 .Xaddr ,Yaddr-1  ,Ls ide  ,Kond  ) 1  node  t  lower 

1  IF  Nd1<>Nd2  AND  Ndl >0  AND  Nd2>0  THEN  1  Evaluate  along  side  nodes 
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!  LET  Tmp I -Oiel< Pixl( Xcnt  ,Ycnt  +  1  > > t  Upper  dielectric  pixel 

!  LET  Tmp2-Diel ( Pix 1 ( Xcnt .Ycnt )> l  Lower  dielectric  pixel 

»  LET  Vt 1 -Hpi v<  Ndl  ) 

!  LET  Ut2-Hpiv( Nd2 ) 

'  IF  Tmp1O0  AND  Tmp2<>0  THEN  Tmp-(  Ut  1  *Tmp  1  +Ut2*Tmp2  )/ (  Tmp  1  +Tmp2 

!  LET  Pot nt ( Xcnt  , Ycnt  + 1  )-< Ut 1 +Tmp  )/2 1  Pixl  volts 

t  LET  Potnt( Xcnt  , Ycnt  )=( Ut2  +  Tmp >/2 

i  IF  Tmpl <>0  AND  Tmp2<>0  THEN 

1  LET  Dsp lc( Xcnt  , Ycnt  + 1  )*( Ut 1 -Ut 2  )/( 1 /Tmp 1 +1 /Tmp2 > 

i  END  IF 

1  LET  Dsplc( Xcnt  ,Ycnt  )-Dsplc( Xcnt ,Ycnt  +  1  ) 1  Displacement  mag. 

i  END  IF 

i  NEXT  Nd 

i  NEXT  Ycnt 

i  END  IF 

1  MAT  Potnt0*  Potnt»(Resp)  1  Normalizing  to  1  volt  across  sample 

1  MAT  Dsplc-  Dsplc*(Resp)  1  &  sum  of  displacements  along  row-diei 

1  LET  Ndt-1  !  sign  provider  for  following  loop 

!  FOR  Xcnt-I  TO  Lside 

!  LET  Resp2»Resp2+Nd1 »( Dsplc( Xcnt .Boxes )-Dsplc( Xcnt ,Boxes+ 1  >  ) 

!  LET  Ndl— Ndl 

!  NEXT  Xcnt 

!  LET  Resp2-Resp2/2  !  dielectric  response  perp  to  E 

•***>  Should  be  end  of  calculations,  printouts  follow 
!***>  Printout  of  the  dielectric  pixel  array 
IF  NOT  (Sprss)  ANO  Lside0  3  THEN 

PRINT  "DIELECTRIC  PIXEL  ARRAY.  2-dimensional  jLsides "by" sLside 
FOR  Xcnt-Lside  TO  1  STEP  -1 
FOR  Ycnt-1  TO  Lside 

PRINT  USING  "DDOO.O  ,f ” sPROUND(Diel(Pixl( Xcnt, Yen t)),-1  > 

NEXT  Ycnt 
PRINT 
NEXT  Xcnt 
END  IF 

>***>  Printout  of  the  hopper  array 

i  *  *  #  *  >  Find  the  senes<->paral  lei  factor 
LET  Tmp 1 -FNWnr<  Di e 1 ( *  >  ,Frpx( * ) ,Resp ,Tmp ,9  > 

PRINT 

PRINT  “Composite  Dielectric  Response  Tensor  Components:" 

PRINT  “  principal-" sResp 

PRINT  “  &  ser ies<->paral lei  factor  -"s 

PRINT  PROUNDUmpI  , -3 >»“<+/-" sPROUNDt Tmp , -4 ) s "X  iteration  error)" 

PRINT 

t  IF  Sprss-0  THEN 

1  PRINT  "PIXEL  UOLTAGES  2-dimensional sLsides "by" sLside 

•  Tmp-FNMatprnt ( Potnt (*) ,-Lside > 

1  PRINT  "PIXEL  DISPLACEMENT  FIELD  MAGNITIUDES  , " sLs 1 de f “by " 1 Ls ide 

1  Tmp-FNMatprnt ( Dsplc< *) ,-Ls 1 de ) 

1  ENO  IF 

!#*»>  NOTE:  Tranpose  used  then  it  is  an  additional  cycle  to  Rptr 
Graf(Grpt  ,1  )-Rptr  +  Tmp/ 1 00 
Graf(Grpt  ,2  )-Frpx( 1  ) 

Graf(Grpt ,3)-Resp 
Graf ( Grpt  ,4  )«Tmp 1 
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LET  Grpt«Grpt+1  !  Increment  storage  counter 

i NEXT  Trans 
NEXT  Rptr 

!***>  output  repeat  calculations 

LET  Bhdr$-"< “&UAL*( Lside >&Hx“&VAL*( Lside )&"  )* 

IF  Tls-1  THEN  LET  Bhdr$-Bhdr*& “  elmnts" 

IF  Tls-2  THEN  LET  Bhdr*-Bhdr*&“ / ( 2x2s > " 

IF  Kond-1  THEN  LET  Bhdr$-Bhdr*S"  InslBC" 

IF  Kond-2  THEN  LET  Bhdr*-Bhdr*S"  PrdcBC" 

LET  Bhdr$-BhdrS& *  Sparse"  !  solution  by  sparse  methods 
IF  Qdrnt-2  THEN  LET  Bhdr$-Bhdr$&“  4fold" 

IF  Ptrn-0  THEN  Bhdr$-Bhdr*&“  intrnl," 

IF  Ptrn-1  THEN  8hdr$-8hdr$&"  “&FlnS 
IF  Ptrn-2  THEN  Bhdr$-Bhdr*i"  USER," 

IF  Ptrn-3  THEN  BhdrS«8hdr$&“  RANDOM,* 

IF  Ptrn-4  THEN  Bhdr$-Bhdr$8.“  SLANT," 

IF  Ptrn-5  THEN  Bhdr$-Bhdr*& “  ELLIPSE," 

IF  Ptrn-6  THEN  Bhdr$-Bhdr$&"  STRAT  , " 

IF  Ptrn-7  THEN  8hdr*-8hdr$&"  BOXES," 

LET  Occp«LEN(Bhdr*) 

LET  Bhdr*I 1 +0ccp 1-RPT$( “  ",80-Occp)  -1  pad  with  blanks 

LET  BhdrSI  601“"  "8,DATE$(  TIMEDATE  >&"  .  "&TIME*(  TIMEDATE  ) 

LET  Dhdr*-Bhdr* 

IF  Rep-1  THEN  PRINT  "  for  the  case  abbreviated  . . " 

IF  Rep-1  THEN  PRINT  Bhdr$ 

IF  Rep> 1  THEN 

PRINT  “  Summary  of “ iRep s "repeat  variations:  (as  programmed)" 

FOR  Rptr-1  TO  Rep 

PRINT  “  Case" i< (Rptr-1  )  DIV  Pose >+1 ;">"  ,PR0UND( Graf < Rptr  ,  1  >  ,-3  ) , 
PRINT  PR0UND< Graf ( Rptr  ,2  >  .-3  )  ,PR0UND<  Graf <  Rptr  ,3  )  ,-3  > . 

PRINT  PROUND( Graf ( Rptr  ,4  )  ,~3  > 

NEXT  Rptr 

OISP  “  Save  repeat  info  (array  form  S I ZE  < Gra f  .  1  )  i "x " i SI ZE( Graf  ,2 > 

INPUT  “)?  0)  No  1)  Definitely" ,Nd1 
IF  Ndl -1  THEN 

DISP  “  Enclose  (in  """“s)  new  file  name  to  send  info  vectors  to7" 
INPUT  "  (null-use  old  file)*, Fin* 

IF  P0S( Fin*,":"  )-0  THEN  Fln*-Fln*&Msd* 

INPUT  "  Title,  (up  to  80  characters  )“  ,AhdrS 

LET  Andril 1 +LEN( Ahdr*  )  ]-RPT*( "  "  ,S0-LEN( AhdrS >  ) 1  pad  with  blanks 
DISP  *  File  named  “““jFIn*i""B  ( [ " ;LEN< Fin*  ) ; " 3  characters  )’ i 
DISP  "  to  contain  repeat  info" 

PRINT  “  File  “ “ “ jFln*j “ “ “ ’ s  user  and  description  headers  are  “ i 
PRINT  "(2  lines):" 

PRINT  Ahdr* 

PRINT  Bhdr$ 

IF  Fin*-""  THEN 

INPUT  “  Enter  the  filename  to  be  created7  nul 1-stop" , Fin* 

IF  Fin*-""  THEN  STOP 

OISP  "  Enter  f  l  le”  ”  ’  ;Fln*» *  *  * ’ s  storage  size  limit  in  bytes  <*“ 
DISP  UALSt  256+8*Rep*Pose»Kwd ) i ” )“ i 
INPUT  "7"  ,Nd1 

IF  Ndl < 1048  THEN  Ndl-1048  •  1  kiloBYTE  run 

CREATE  FlnS  ,Ndl 
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6930  ELSE 

6940  IF  POSlFln*  : “ >-0  THEN  FlnS-FlnSSMsd* 

6950  END  IF 

6560  ASSIGN  Slnfostr  TO  Fin* iFORMAT  OFF 

6970  OUTPUT  @Inf ostri Ahdr* ,8hdr$ ,Rep ,Kwd ,6raf < * >  .END 

6980  AS5I6N  SNodstr  TD  » 

6950  END  IF 
7000  ENO  IF 

70(0  (#**>  ptxel  file  output  choice 

7020  LET  Nd ! -0 

7030  IF  Ptrnoi  THEN  INPUT  “  Save  last  pixel  grid?  0  )No  1  )Yes"  ,Ndl 
7040  IF  Ndl-i  THEN 


70s0  DI5P  *  Enclose  (in  *“*“s)  new  file  name  to  send  pixel  pattern  to"'*; 
7060  INPUT  “  inuii-use  old  fiie>“,Fln* 

7070  IF  P0S(Fln*  >-0  THEN  F  InS-F  lnSS-Msd* 

7060  INPUT  '  Title  (up  to  80  characters  )  if  null  then  default  label "  ,Chdri 
7090  LET  0hdr*-8hdr» 

7100  OISP  *  File  named  * " “ sFlnSi " * “  ( C " »LEN< Flni  ) ; " ]  characters)"! 

71(0  OISP  "  contains  the  pixel  grid" 

7(20  PRINT  *  Fiie  “  “ “ ;F In* ; “ “  “ 1 s  header  is  “ 

7130  PRINT  Hdr* 

7140  IF  FlnSO"  ”  THEN 

7i50  DI5P  "  Give  f lie'"”  ;FIn*i *,*’s  max  capacity  limit  in  bytes*; 

7 1 60  DI5P  “?  (” “ ;UAL*< (28+SHIFT<Px_tot  ,-(  ));“>* ; 

7i 70  INPUT  “  *  ,Nd1 

7 i 80  IF  Ndl <256  THEN  Ndl-256 

7i90  CREATE  Fln*,Nd( 

7200  ElSE 

72(0  INPUT  *  Enter  the  existing  f l lename? “  ,Fln* 

7220  IF  P0S( Fin*  , “ : “  )-0  THEN  Fln*-Fln*&Msd* 

7230  END  IF 

7240  ASSIGN  @Pixstr  TO  Fin*; FORMAT  OFF 
7250  OUTPUT  @Plxstr;Hdr*  .Lside ,Pixl( * ) ,ENO 

7260  ASSIGN  @Pixstr  TO  * 


7270  END  IF 

7280  !***>  Interaction  file  output  choice 

7250  PRINT  RPT*( *  *  ,25  );"... elapsed” ; PROUNDi TIMEDATE-Star i  i  ) ; 

7300  PRINT  “sec  for  completion  at  “ ; TIME*( TIMEDATE  ) 

73(0  PRINT  "  MEMORY  IS* ;UAL( SYSTEM*! "AVAILABLE  MEMORY”  )  >/8 » * t reals ) ” 

7320  LET  Lxtnt-Lside  *  update  COM  /Pixei/  le  Pixl<*>  size 

7330  END 

7340  •  n  H  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it 

7350  DEF  FNNi( INTEGER  Xn  ,  Yn  ,Lszn  ,Lup  ) 

7360  !  Returns  the  node  number  for  a  square  Pixel  grid  network 

7370  !  of  capacitors  for  case  of  insulated  sides. 

7360  t  (Xn,Yn)  2D  coordinates  of  Pixel  leading  to  the  nearest  node 
7350  !  If  Up-1  then  Pixl  above  node,  else  Pixl  below  node 

7400  INTEGER  Xhf  ,rsw  ,Ndb  ,Lyr 
74 i0  LET  Lup-6IT( Lup  ,0 ) 

7420  LET  Lyr-SHIFT( Lszn ,  I  )  1  in  essence  divides  by  2 

7430  IF  Xn>0  AND  xnOLszn  AnO  Yn>0  AND  Yn<-Lszn  THEN 

7440  LET  Ndb-1 

7450  IF  Kn-1  AND  Lup  THEN  LET  Ndb-0 

7460  IF  Xn-Lszn  AND  NOT  (Lup)  THEN  LET  Ndb-5HIFT( Lszn»Lszn  ,  (  )-Lszn+2 
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7470 

7480 

7490 

7500 

75  ra 

7520 

7530 

7540 

7550 

7550 

7570 

7580 

7590 

7600 

7610 

7620 

7630 

7640 

7650 

7660 

7670 

7680 

7690 

7700 

7710 

7720 

7730 

7740 

7750 

7760 

7770 

7760 

7790 

7600 

7810 

7820 

7650 

7840 

7650 

7860 

7670 

7880 

7890 

7900 

7910 

7920 

7930 

7940 

7950 

7960 

7970 

7980 

7990 

8000 


IF  Ndb- I  THEN  *  executes  if  node  not  yet  assigned 

LET  Xhf -SHIFT! xn- 1 -Lup  ,  1  )  I  effectively  Dlv  2 
LET  Xsw-( BIT< Xn  ,0  )  EXOR  Lup  H  0-even  1-odd  in  Xn 
LET  Ndb-Xhf *< Lszn-1  >+5HIFT< Yn+Xsw  ,  I  H  basic  node  count 
IF  NOT  <  Xsw )  THEN  LET  Ndb-Ndb+Lyr 

IF  (Yn-t>  AND  NOT  (Xsw)  THEN  LET  NdP-Ndb+Lyr-Lup*<  Lszn-U 
IF  Yn-Lszn  AND  NOT  (Xsu)  THEN  LET  Ndb-Ndb+Lyi — ) -Lup* i Lszn- 1  ) 

END  IF 
ELSE 

LET  Ndb-- I 
END  IF 
RETURN  Ndb 
FNENO 

1  ][  U  it  H  ][  ][  H  H  H  H  H  K  ][  3C  H  H  JC  H  ][  H  K  H  If  1C 
OEF  FNNp( INTEGER  Xn  ,Yn  ,Lszn  ,Lup  ) 

!  Returns  the  node  number  for  a  square  Pixel  grid  network 
i  of  capacitors  for  case  of  periodic  or  wrap  around  sides 
i  (Xn.YnJ  20  coordinates  of  Pixel  leading  to  the  nearest  node 
1  If  up-l  then  Pixi  above  node,  else  Pixl  below  node 
INTEGER  Ynn.Ysw.Ndb.Lyr 
LET  Lup-BIT(Lup  ,0  ) 

LET  Lyr-SHIFT< Lszn  ,  1  )  i  in  essence  divides  by  2 

LET  Ynn- i +( (Yn- 1  )  MOD  Lszn)  !  wrap  around  in  Y  coordinate 

IF  Ynn<  t  THEN  Ynn-Y'nn+Lszn 
IF  Xn>0  AND  Xn<-Lszn  THEN 
LET  Ndb- I 

IF  Xn-1  AND  Lup  THEN  LET  Ndb-0 

IF  Xn-uszn  AND  NOT  (Lup>  THEN  LET  Ndb-SHIf T( Lszn*Lszn  ,  1 )-Lyr+l 
IF  Ndb- 1  THEN  !  executes  if  node  not  yet  assigned 

LET  X'sw-(  BIT(  Xn  ,0  )  EXOR  Lup  H  0-even  I -odd  m  Xn 
LET  Ndb«< Xn- l-Lup  )*Lyr+5HIFT< Ynn+Xsw , I ) f  basic  node  count 

IF  Ynn- t  AND  nOT  (Xsw)  THEN  LET  Ndb-Ndb+Lyr 
END  IF 
ElSE 

LET  Ndt>-- 1 
END  IF 
RETURN  Ndb 
FNEND 

•  1C  1C  If  K  If  If  H  If  H  H  If  H  If  ]I  1C  If  If  1C  H  If  H  If  If  1C 
SUB  C'vndi<  INTEGER  Nodi  ,Lszi  .Sctri  ,Xout  i  ,Youti  ) 

1  Converts  a  node  number  of  layering  scheme  into  the  (x,y> 

1  coordinates  of  neighooring  nodes 
!  IN  Nodi  -  node  number 

’  IN  Lszi  -  Pixl  extent  along  either  X  or  Y 

!  IN  Sctri  -  selection  adjacent  Pixl  neighbor  to  node 

•  (0-  X  lower, Y  iowerj  1-  X  lower, Y  higheri 

'  2-  X  higher,  Y  lowen  3-  X  higher,  Y  higher) 

1  OuT  xouti  -  X  coordinate  address  outcome 

i  OUT  Youti  -  Y  coordinate  address  outcome 

1  internal  vanaoles: 

1  Lyri  -  a  counter  for  number  of  layers 
1  Swi  -  an  even  odd  switch 
1  Ovri  -  overfill  counter 

1  Nodmaxi  -  maximum  node  for  given  Pixel  grid  size,  Lszi 
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8010  i  Hfszi  -  half  of  Pixel  grid  size,  Lszi  (then  Lszi  must  be  even) 

80Z0  INTEGER  Lyri  ,5wi  ,Gvri  .Nodmaxi ,Hf szi 

8038  IF  SIT(Lszi,0)  THEN  PRINT  *  warning,  odd  Pixel  grid  extent" 

8040  LET  Nodmaxi-SMIFT(Lszi*L3zi , I ;-Lszi+2 

8050  LET  Hf szi-SHIFT < Lszi  .  1  )  •  effectively  DIV  2  operation 

8060  SELECT  Nodi 

8070  CASE  <1 

8080  LET  Xouti“0 

8080  LET  You t i*0 

8100  CASE  >-Nodmaxi 

8110  LET  Xout i*Lsz i+ 1 

8 i 20  LET  Youti"Xouti 

8130  CASE  ELSE 

8140  LET  Lyn»  i *k  ( Nodi- 1  )  DIU  (Lszi-I  l  )  1  for  Biiayer 

8(50  LET  Ovri* t +( ( Nodi- I  )  MOD  (Lszi-!))  !  for  4  Nodes  in  Bilayer 

8(60  LET  Swi ”( Uvr i >Hf sz l  )  !  0"iower  l*upper  in  Bliayer 

8170  IF  Swl  THEN  lET  0vri“0vrl-Hf szi  *  adjust  for  t  Nodes  in  upper 

8180  LET  Xout i”SHIFT< Lyr i  ,- 1  )-) +Swi+BIT( Sctri  ,  1  ) 

8150  LET  Yout  i“SHIFT(  Ovn  ,-1  )-1 +Swi+BIT(  Sc tr l  ,0  ) 

8200  END  SELECT 

8210  SUBEND 

8220  f  ][  It  It  It  It  It  K  It  It  It  It  H  It  It  It  It  It  It  IE  it  IE  IE  It  IE 
8230  SUB  Cvndp<  INTEGER  Nodp  , Lszp , Sctrp , Xoutp ,Youtp ) 

8240  !  Converts  a  node  number  of  layering  scheme  into  the  (x.y) 

8250  1  coordinates  of  neighboring  nodes 

8260  >  IN  Nodp  “  node  number 

8270  *  IN  Lszp  -  Pixi  extent  along  either  X  or  Y 

8280  i  IN  Sctrp  -  selection  adjacent  Pjxi  neighbor  to  node 

8280  *  ( 0»  x  lower  ,Y  ioweri  l*  X  lower, Y  higher! 

8300  !  2-  X  higher ,  Y  lower f  3-  X  higher,  Y  higher) 

8310  '  OUT  Xoutp  «  X  coordinate  address  outcome 

8320  t  OUT  Youtp  »  Y  coordinate  address  outcome 

8330  i  internal  variables: 

8340  1  Qvrp  *  overfill  counter 

8350  !  Nodmaxp  -  maximum  node  for  given  Pixei  grio  size,  Lszp 

8360  *  Hfszp  »  half  of  Pixel  grid  size,  Lszp  (then  Lszp  must  be  even) 

8370  INTEGER  Ovrp  .Nodmaxp  .Hfszp 

8380  IF  BIT(Lszp,0)  THEN  PRINT  “  warning,  odd  Pixel  grid  extent" 

8390  LET  Hfszp-SHIFTtLszp  ,  I  )  *  effectively  DIV  2  operation 

8400  LET  Nodmaxp*SHIFT( Lszp*Lszp , I  )-Hf szp+ I 

84(0  SELECT  Nodp 

8420  CASE  < I 

8430  LET  Xoutp“0 

8440  LET  Youtp“0 

3450  CA5E  Y“Nodmaxp 

8460  LET  Xoutp-Lszp-H 

8470  LET  Youtp“Aoutp 

8480  CASE  ELSE 

8490  LET  xoutp" t +<( Nodp- I  )  DIV  Hfszp)  1  for  slab  of  X 

6500  LET  Gvrp“ I +( < Nodp- !  )  MOD  Hfszp)  •  for  t  Nodes  within  slab 

85(0  LET  Youtp-SHIFT( Ovrp  ,- I  )+BIT( Sc trp  ,0  )-B IT( Xoutp  ,0  ) 

8520  LET  Xoutp,*Xoutp+8IT<  Sctrp  ,  I  ) 

8530  IF  YoutpYLszp  THEN  Youtp-Youtp-Lszp 

8540  END  SELECT 
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8550 
8560 
8570 
8580 
8580 
8600 
86  10 
8620 
8630 
8640 
8650 
8660 
8670 
8680 
8690 
8700 
8710 
8720 
8730 
8740 
8750 
8760 
8770 
8780 
8790 
8800 
88(0 
8820 
8830 
8840 
8850 
8860 
8870 
8880 
8890 
8900 
89(0 
8920 
8930 
8940 
8950 
8960 
8970 
8980 
8990 
9000 
90(0 
9020 
9030 
9040 
9050 
9060 
9070 
rara 


SUBEND 

i  JC  It  3t  It  It  It  3t  3C  3t  It  3t  H  3C  H  H  3t  It  3C  3t  3C  3C  JC  3f  3E 
DEF  FNMatprnt < Matrxt *),  INTEGER  Ordr  ) 

t*#»>  printout  of  an  array  sized  Ordr  x  Ordr  &  @  element  Spcs  wide 
COrl  /Pass/'  Relay 
INTEGER  Sbcx ,5bcy ,Sbcyy  ,Sbc  ,Typ 

LET  Sbc-SGNt Ordr  )  i  Reverse  printout  of  rows  indicator 

LET  Pslct-0  i  Pslct  is  the  biggest  element  magnitude 

FOR  Sbcy- I  TO  ABStOrdr) 

FOR  Sbcx- 1  TO  ABStOrdr) 

IF  Pslct  <A6S<  Matrxt  Sbcx  , Sbcx  )  )  THEN  LET  Pslct-ABSt Matrxt  Sbcx  , Sbcx  ) ) 
NEXT  Sbcx 
NEXT  Sbcy 

IF  Psict>0  THEN  LET  Typ-2-INTt LGT< Ps let  )  ) 

IF  Typ-4  OR  Typ-5  THEN  Typ-3 
IF  Typ--I  THEN  Typ-0 
FOR  Sbcy-I  TO  ABStOrdr) 

LET  Sbcyy-5bc*Sbcy+<A6S<0rdr)+l  )*<  l-5bc)/2 
FOR  5bcx-(  TO  ABStOrdr) 

SELECT  Typ 
CASE  -0 

PRINT  USING  “DDODD ,i“ ;PROUNOtMatrxtSbcx  .Sbcyy  )  ,0) 

CASE  -I 

PRINT  USING  ”000 . 0  ,#* i PROUNDt Matrx( Sbcx, Sbcyy  ),-l  ) 

CASE  -2 

PRINT  USING  “00.00,*“ j PROUNDt  Matrxt Sbcx  .Sbcyy  )  ,-2  ) 

CASE  -3 

PRINT  USING  “0.000  ,1’’  ?PROUNDt  Matrxt  Sbcx  ,  Sbcyy  )  ,-3  ) 

CASE  ELSE 

LET  Typ-56N< Typ  ) * 99 
END  SELECT 
NEXT  Sbcx 

IF  ABS<  Typ  )<>99  THEN  PRINT 
NEXT  Sbcy 

IF  Typ-99  THEN  PRINT  “  ...  array  too  small  to  format" 

IF  Typ--99  THEN  PRINT  *  ...  array  too  big  to  format” 

RETURN  t 
FNENO 

<  3 1  3t  31  It  31  3t  It  3 1  3t  3t  31  3t  3t  3t  3t  3t  It  31  31  It  It  3t  3C  3t 
DEF  FNWnrt Oiel (*)  ,Frpxt •)  ,0in  .Pcterr , INTEGER  Nth) 

Object  of  function  to  find  alf  satisfying  l-SUM  oft  vol*diel*alf  > 
'***>  where  vol-fract  lonal  volumes 

diel«t species  permittivity  )/( composite  permittivity) 

!••*>  alf-  constant  exponential  between  -I  &  +1  < ser ies<->paral lei  ) 

COft  /Pass/  Relay 
INTEGER  Spk.Arnd 

IF  Oin-0  OR  MAXtFrpxt *  )  >>-l  THEN 
LET  Alfl-0 
LET  Pcterr-t .E-99 

PRINT  “  ..  series  <->  parallel  factor  at  or  beyond  limits” 

ELSE 

LET  Gsum-0 
LET  6dev-0 
FOR  Snk-I  TO  Nth 
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9100 
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9600 

96  1 0 

9620 


LET  brat io*0 

IF  QleiiSpk>>0  ANO  Frpx(Spk)>0  THEN 
LET  Sratio-L06(Diei< Spk >/Uln ) 

LET  Gsum-Gsum+Frpxl Spk  )*Grat io 
LET  Gde v-bdev+Frpx ( Spk  )*bratio*Gratio 
END  IF 
NEXT  Spk 

LET  A1 f 1 --2 »bsum/'bdev 

LET  bsum-bsum+bdev  i  Initiai  guess  of  sun 

LET  Try-1  l  Optimize  sign  in  iterations 

LET  Aif0-0 
LET  Arnd-1 

WHILE  AB5(Alf!-Alf0)>. 0000000001 
LET  bdev-0 
LET  6sum*0 
FOR  Spk- I  TO  Nth 

LET  Gra 1 10-D lei < Spk  )/Din  1  relative  permittivity  to  composite 

LET  6wk-0  !  a  term  @  species  to  sum 

IF  Gratio<>0  AND  Frpx(Spk)>0  THEN 
LET  Gwk-Frpx(  Spk  )*6rat  io"Al  f  | 

LET  Gsum-Gsum+Gwk  ' -  sum  function  (as  described) 

LET  bdev-6dev+LGG( Grat  10  )*Guik  1  1st  derivative 
ENO  IF 
NEXT  Spk 
LET  Pcterr-bsum 
LET  Gsum-Al  f  I  *  (  Gsum- 1  )/'6dev 
LET  Gdev-Al f I *A1 f 1 -Try*6sum 
IF  6dev<0  THEN  Gdev-Gdev+2*Try*6sum 
LET  Alf2-SGN(Alf I  >*SQR(6dev) 

IF  Arnd  MOD  2-0  AND  ABS< A1 f2-Al f I  )>ABS< A1 f 1 -A1 f0  )  THEN  Try— Try 
LET  Aif0*»Alft  !  bumping  iterations 

LET  Ai f  t-Al f 2 

LET  Arnd-Arnd+I  !  incrementor 

END  WHILE 

IF  AbslAif  K. 000000 1  )  THEN 
LET  Pcterr- 1 00* < Pc  ten — I) 

ELSE 

LET  Pc terr- 1 00»< Pcterr  '< l/AIfl  )-l  ) 

END  IF 
END  IF 
RETURN  Alft 
FNEND 

i  H  K  It  ] i  H  j (  1 C  3 C  H  3 [  3 1  3 i  3 i  31  3 i  3 i  3i  3 i  3 i  3 [  3 i  3 i  3 i  31 
SUB  Pixi2d  tlit 


'  This  program  is  for  patterning  a  dielectric  bond  grid  (rectangular 
1  type)  with  a  binary  mixture.  The  interface  is  approximately  flat 
i  at  a  tilted  sngie.  (Say  like  a  tilted  giass  of  water) 

i  + - + 

'  1 . .  1 

1  1  ..  1 

'  i  . .  1 

,  + - + 
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5630 
9640 
5650 
5660 
5670 
5650 
9690 
9700 
57 10 
9720 
9730 
9740 
9750 
9760 
9770 
9780 
5790 
9800 
5810 
9820 
9830 
5840 
9850 
5860 
5870 
9880 
9890 
9900 
5910 
9920 
9930 
9940 


DE6 

COM  /Pass'  Relay 

COM  /Pixel/  Chdr$C 80 1  .DhdrSC 80 ], INTEGER  Lbonda  ,Pix 1 ( I : 1 80 , 1 : 1 80 ) 
INTEGER  Accum  ,0pt  , Ip  ,  Jp 
REDIM  Pixlt  t  tLbonds  ,  I  :Lbonda  ) 

PRINT 

PRINT  *  Makes  an  bond  grid  pattern  with  a  binary  mixture  whose 
PRINT  "interface  levei  sits  at  a  selected  volume  fraction  &  tilt.* 
INPUT  "Request  the  desired  volume  percentage  of  tne  1st  component "  ,Cmp 
InPuT  request  the  desired  tiit  angie  < in  degrees r;  ,Degt 1 t 
PRINT  "The  requested  volume  ls"sCmp;"7.  &  tilt  angie  l  s  *  ;  Degt  i  t ;  *  deg" 
PRINT  “for  a  square  grid  of  bonds  of  size* iLbonds i “X* ;LDonds \ " . " 

OISP  “Decide  fate  of  straddlers  “ ; 

INPUT  " 1  >  pixel  pop/ volume  demand  2)  individual  volume  sway*  ,0ptn 
LET  Cmp«Cmp/100 
LET  Pmc-l-Cmp 

i  **»»••**»•*»••*•*•*•»****»»••••• 


•  From  the  information  of  the  two  components  &  the  angie  of  tiit 
1  then  interfacial  level  can  be  drawn  across  the  grid  window. 

(»•*•*•*»••**«**»***•»»»*»•«••»« 


LET  Sn«A85< SIN( Degt 1 1  )  )  ! 

LET  Cs“AB5< COSt Degt 1 t J >  ! 

LET  Crit“MIN( Sn  ,Cs  )/MAX<  Sn  ,Cs  >/2 1 
Dgmin-MIN< 5n  ,Cs >/< 5n+Cs >  ! 

LET  Ogsn-Sn/ ( 5n+Cs >  ! 

LET  Ogcs-t-Dgsn  t 

SELECT  Cmp  ! 

CASE  <>Crit 


LET  Cpt-SGR< 2*Cmp*0gsn*Dgcs  i 
CASE  >-( 1 -Cri t  ) 

LET  Cpt»i  — SCjR t  2*Pmc*Dgsn*Dgcs  ) 
CASE  ElSE 


5me 

Cosine 

Cricticai  volume  transition 
A  transition  on  the  diagonal 
Diagonal  threshold  at  Cnt  vol 
&  its  complement 
y-intercept  for  different  cases 


9950 
9960 
5570 
9980 
5990  • 

1 0000  I 
10010  I 
[0020  I 

10030  I 
(0040  I 
10050  I 
(0060 
(0070 
1 0080 
10050 
(0(00 
(0M0 
(0(20 
10130 
10140 

1  0150 
1  0160 


lET  Cpt m( Cmp- .5  )/ u +2*Cr i t )+  .5 

END  SELECT  i  Done  computing  diagonal  intercept 

'PRINT  “The  level  has  slope" iPROUNDl Sn/Cs ,~2 );  “&  a  diag-intercept  of"; 
'PRINT  PROUNDlCpt ,-3  >;" . * 


The  array  wiii  be  filled  according  to  whether  the  element 
happens  to  be  on  one  or  tne  other  side  of  tne  tiit  iine. 
If  an  element  is  crossed  through  Dy  tne  tiit  line  then 
the  element  will  be  assigned  to  the  componet  possessing 
the  greatest  volume  deficet  as  the  pixels  are  filled. 


i  Initialize  keep  tabs  of  excess 
1  Initialize  totals  counter 


LET  vbias»0 
LET  Accum«0 
FOR  Ip«(  TO  Loonds 
FOR  Jp-t  TO  Lbonds 

0iag"Jp*0gsn+Ip*0gcs- .5  <  Diag  intercept  ©pixel 

Di f "Diag-Cpt^Lbonds  '  Distance  offset  on  diag 

REM  PRINT  PROUnO( 0  lag  ,-2  ) ! PR0UND< D l f  , -2  ) ;  1  Matrix  calc  out 


SELECT  Olf 
CASE  <■- .5 


Pixl( Ip  ,jp 
CASE  /-.5 


1  Assign  1st  component 
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10170 
10180 
10190 
10200 
10210 
1 0Z20 
i  0230 
10240 
10250 
10260 
10270 
1 0280 
10290 
10300 
10310 
i  0320 
i  0330 
10340 
10350 
10360 
(0370 


Fix  1 ( Ip  ,  Jp  )-2 
CASE  ELSE 
Dpf-.S-Dif 
SELECT  Dpt 
CASE  <Dgmin 

Upx 1-Dp  f *Dp  f /Ogsn/Dgcs/2 
CASE  X  I -Dgmin  ) 

Upx 1- I -( <  t  -Dp  f  )  '  2  )/'Dgsn/Ogcs/2 
CASE  ELSE 

Upx 1 -( Dp f - . 5  )* (  I +2«Cr i t  )  + .  5 
END  SELECT 

REM  PRINT  PROUNOt Opxi  ,-2  )  i 
LET  Vbias-Ubias+Upxl 
IF  Optn-2  THEN  Ubias-Vpxl 
IF  Vb i a s > . 5  THEN 
Fixlt Ip  ,Jp  )  - 1 
vbias-Ubias-  I 
ELSE 

Pixli  Ip  ,Jp  )-2 
END  IF 
EnD  SELECT 


•  Assign  2nd  component 
1  Borderline  pixel  case 
1  Diagonal  witnin  pixel 


1  Uoi  fra c  u/ittun  a  pixel 
•  For  matrix  vanoies  printout 
i  Accum  of  excess  component  1 
I  Option  of  straight foward  bias 


10380  Accum-Accum+2-PixK Ip  , jp )  1  Accumulation  of  1st  component 

10390  NEaT  Jp 
10400  N£XT  Ip 
10410  SUBENO 


10420  !  ][  K  It  It  II  II  31  31  31  31  51  3i  31  H  31  3i  H  K  3i  31  31  31  31  3i 
10430  SUB  Pixi2d_rand 

10440  !***>  Subprogram  to  create  a  ZD  pixel  array  of  random  distribution 
10450  !  "&%  )  l  RANDOM  2D  PIXEL  GRID  it ! »X >< V $* 

10450  COM  /Pass/  Relay 

10470  COM  /Pixel/  Chdr*C 80 ] ,0hdr*[ 80 J . INTEGER  Lpix  ,Pixl< t : I  80 , 1 : I  80 > 

10480  INTEGER  Xp ,Yp ,Xq ,Yq .Fill ,Frdm ,Sqrs ,Pxs .Pixtmp  .When 
10490  REDIM  Pixl< I :Lpix  ,  t JLpix  ) 

10500  PRINT  ‘  Enjoy  creating  a  randomized  2D  pixel  grid  uinose  pixei  elements' 
10510  PRINT  are  iabeiled 

10520  INPUT  “Random  seed7  (integer  or  neg  if  to  be  via  timer!"  .Vchk 
10530  'LET  Uch'k -Relay 

10540  IF  Uchk-0  THEN  RANDOMIZE  INT ( TIMEDATE  MOD  32757; 

10550  IF  Ucnk>0  THEN  RANDOMIZE  INT(Vchk) 

10560  LET  Fiii-0 

10570  LET  5qrs“Lpix*Lpix 

10580  LET  Frdm-0 

(0590  LET  When-0 

10600  LET  Pxs-0 

10610  LET  Yp-I 

10620  LET  Uchk-0 

10630  'let  Rq- !00*Relay/64  '  if  to  program  using  the  COM  Relay 

10640  FOR  Fill- I  TO  Sqrs 
10550  IF  Pxs<1  THEN 
*0660  LET  When-When+ 1 

10670  OISP  ‘Filling  with  component  I “ ;UAL$< When M *  1  ,  give  volume  "  s 

10680  INPUT  ‘percent ‘  ,Rq 

10690  IF  Rq<0  THEN 

10700  LET  Pxs-Sqrs-Fill+I 
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ELSE 

LET  Pxs-INT< Rq»Sqrs» .01  +  . 5  ) 

IF  Pxs/Sqrs-F 1 1 1  + 1  THEN  LET  Pxs“Sqrs-Fi  1 1+  1 
ENO  IF 

PRINT  '  Component  C  1 5UAL$l When  )  i  ‘  1  is  assigned" iPxs; ‘pixels' 

Enu  IF 

LET  Xp*Fill-Frdm 
IF  Xp>Lp  1  x  THEN 
LET  Xp«Xp-Lpix 
LET  Vp«Vp+1 
LET  Frdm-Frdm+Lpix 
END  IF 

IF  Pxs/0  THEN  LET  P ix i ( Xp  ,Yp  )*Wnen 
Pxs“Px3- i 
NEXT  Fill 
PRINT 

(***■>  iotto-ing  or  random  nixing 
FOR  Fili*l  TO  5qrs 

LET  Frdm“INT<  I +RND»Sqrs  ) 

IF  FillOFrdn  AND  Frdm<“5qrs  THEN 
LET  Xp- I +( <  F 1 1 1 - 1  )  MOD  Lpix) 

LET  Vp- 1 +< (Fill-1  >  DIO  Lp ix  ) 

LET  Xq« I + ( ( Frdm- I  )  MOD  Lpix) 

LET  Yq-l+< (Frdm- I  )  DID  Lpix) 

LET  P ix tnp*P ix 1 < Yp  ,Xp  ) 

LET  Pixl( Yp ,Xp >«Pixl( Vq  ,Xq  ) 

LET  PixK  Yq  ,Xq  J-Pixtmp 
END  IF 
NEXT  Fill 
SUBEND 

'  It  1C  H  It  H  H  It  1C  1C  ][  H  K  It  1C  It  It  1C  1C  It  1C  IC  It  1C  It 
SuB  Pixl2d_fill 

(*•♦>  Subprogram  to  output  a  pixel  array  hand  filled  by  user 
COM  /Pass/  Relay 

COM  /Pixel/  Chdr$t  80 1 ,0hdr$C 80 1 ,  INTEGER  Lpix ,Pixl( ( : 180 , 1 : 180  ) 

INTE6ER  Xpix ,Ypix .Fill  ,Xp  ,Yp 
PEDIM  Pixl( I :Lpix  ,  I SLpix  ) 

MAT  Pixl-  (0) 

FOR  Xpix*l  TO  Lpix 
FOR  Xpix-t  TO  Lpix 
CLEAR  SCREEN 
PRINT 

FOR  Xp-Lpix  TO  I  STEP  -I 
FOR  Yp- I  TO  Lpix 

PRINT  USING  “0D  iPixl< Xp  ,Yp  ) 

NEXT  rp 
PRINT 
NEXT  Xp 
PRINT 

DISP  "Filling  at  i  “  i  xpi  x \ -  ,  t tp  1  x  r:  1  “ ; 

INPUT  r  specify  species  type  witn  integer  1  .  .  9‘ .Neuipix 
LET  Newpix-Newpix  MOD  10 
LET  P ix  1  ( Xpix  ,Ypix  )-Newpix 
NEXT  Ypix 
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1250  NEXT  Xpix 

1260  CLEAR  SCREEN 

1270  FOR  Xp-Lpix  TO  I  STEP  -I 

1280  FOR  Yp-t  TO  Lpix 

I  2S0  PRINT  USING  DO  ;Pixl( Xp  ,Yp  ) 

1300  NEXT  Yp 

1310  PRINT 

1320  NEXT  Xp 
(330  SU8EN0 

1343  !  H  K  It  H  H  ][  1C  H  1C  It  It  It  H  H  It  It  It  It  It  It  It  It  It  It 
1350  SUB  Pix I2d_el lps 

1360  I  >  >  >  Subprogram  to  fill  a  pixel  window  with  an  ellipse  \  <  <' 

1370  COM  /Pass/  Relay 

1380  COM  /Pixel/  Chdr$[ 80 1  ,Ohdr$t 80 1 , INTEGER  Lwndw  ,Pix 1 ( 1 : I  80 , 1 : I  80  ) 

1390  INTEGER  Opt  ,Xp  ,Yp 

1400  REOIM  Pix l ( I : Lwndw  ,  I : Lwndw  ) 

1410  PRINT  "Choice  of  filling  pixel  grid  with  an  ellipse" 

1420  OISP  “Uolume  percentage  of  C t Is  to  be  filled  by  ellipse"; 

1430  INPUT  “  (if/. 5  then  host  )“  ,Ulm 

I  440  LET  Ulm-Ulm/ 100 

1450  DISP  "Type  I  or  2  to  opt  for 

1460  INPUT  "I)  eccentricty  or  2)  axis  ratio”, Opt 

1470  OEG 

1480  IF  Opt-I  THEN 

1490  INPUT  ”  Eccentricity’”  ,Ecn 

1500  IF  AB5<  Ecn  >> I  THEN  PRINT  "  Too  eccentric” 

1510  LET  Mnjr-SQR<  I -Ecn*Ecn  ) 

1520  ELSE 

1530  INPUT  ”  Ratio  (minor  axis)/(major  axis)?".Mnjr 

1540  IF  A8S< Mnjr  )> I  THEN  Mnjr-I/Mnjr 

1550  LET  Ecn«SQR< t -Mn jr*Mnjr  ) 

1560  END  IF 

1570  INPUT  ”  Angle  of  major  axis  w.r.t.  horizontal  in  degrees?”  ,Mang 
1580  PRINT  "The  ellipse  character l st ics  are:" 

1590  PRINT  “  eccentrici ty*" jPROUNDt Ecn  ,-4  ) ; 

1600  PRINT  ”,  axis  rat io»” ;PR0UND( Mnjr  ,-4  ,  &“ 

1610  IF  Ulm> . 5  THEN  LET  Ulm«.5-Ulm 
1620  LET  Major«SQR( ABS( Ulm  )/Mnjr/PI  ) 

1630  LET  Minor«Mnjr*Ma jor 

1640  PRINT  “  measures  a  sexy " ;PROUND( 2»Minor  ,-4  )  ; 

(650  PRINT  "by " i PROUNDt 2*Ma jor ,-4 ) ; " w .r . t  fraction  of  pixel  window." 

1660  FOR  Xp»l  TO  Lwndw 

1670  LET  Xpos“( ( Xp-Lwndw  DIO  2)~.5>/Lwndw 

1680  IF  A6S(Xpos>>.5  THEN  PRINT  *  Warning.,  touches  window  edge" 

1690  LET  Xcs“Xpos*C05< Mang > 

1700  LET  Xsn"Xpo3*SlN< Mang  ) 

1710  FOR  Yp*l  TO  Lwndw 

1720  LET  Ypos-( < Yp-Lwndw  DID  2)-.5>/Lwndw 

1730  IF  ABS(Ypos)/.S  THEN  PRINT  "  Warning.,  touches  window  edge" 

1740  LET  Xpos"Xcs-Ypos»SIN( Mang > 

1750  LET  Ypos« Ypos*C0S< Mang  i  +  Xsn 

1760  LET  Gf ct-< Ypos/Major  )  2  +  ( Xpos/Minor r 2 

1770  IF  Gfct/I  THEN  Pixl( Xp ,Yp )-2-( l-SGN; Ulm M/2 

1780  IF  6fct'.-l  THEN  PixKXp.Yp  )-!+(!  -SGN(  Ulr  )  )/2 
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1  I  790  NEXT  Yp 
I  1600  NEXT  Xp 
11810  SUBENO 

11820  •  H  It  It  It  H  It  3 C  3t  3t  3 1  3t  3 1  3 1  3t  3t  3t  3t  3t  3t  3t  3t  3t  3t  )t 
11830  SUB  Pixl2d_strat 

11640  i  Pixel  griddmg  of  stratified  layers 

11850  COM  /Pass/  Relay  !  for  sharing  to  subs 

I  1 880  COM  /Pixel/  ChdrSC  80  3  ,0hdr$C  80  3  .INTEGER  Lxtnt  ,Pixl< I s 180 , 1 : I  80  ) 

11870  INTEGER  Strat .Stack  .Fix  .Lent , I s  ,Ns  .Swap 
11880  REOIM  Pixl(  I'.Lxtnt  .1  :Lxtnt  ) 

11890  INPUT  “Layer  sequencing?  random*!  by  user-2 “  .Strat 
11900  IF  Strat <0  THEN  STOP 

11910  LET  Strat -( Strat >  I  )  *  reassigned  0  or  I 

I  1920  IF  NOT  ( Strat  >  THEN 

11930  INPUT  “Random  seed?  (give  0  via  timer  or  give  integer  )" .Seed 
I  1940  IF  Seedx0  THEN  STOP 

11950  LET  Swap*lNT( Seed* . 5  ) 

I  I960  IF  Swap-0  THEN 

11970  RANDOMIZE  INT ( TIMEOATE  MOO  32787) 

I I  960  ELSE 

11990  RANDOMIZE  Swap 

12000  END  IF 

12010  END  IF  !  Strat  test  RANDOMIZE 

12020  DISP  "Orientation  of  layers?  flatwise*-!  spanned-*!  (or  series<->“» 

12030  INPUT  “paral lei  )"  .Stack 

12040  LET  5tack-( Stack>-0 )  !  reassigned  0  or  1 

12050  IF  Strat  THEN 

12060  LET  Is-I 

12070  FOR  Lent- I  TO  Lxtnt 

12080  LET  Ns-0 

12090  DISP  “Assign  to  layer" jLcnt i "a  component  type  [I. .93“? 

12100  INPUT  “  or  null  to  default  to  last”, Ns 

12110  IF  Ns<0  THEN  STOP 

12120  IF  Ns-0  THEN  I  defaults 

1 2150  LET  Ns-Is 

12140  ELSE 

12150  LET  Is-Ns 

12160  END  IF 

I  2  I  70  PRINT  “  Layer " (Lent ! "assigned  to  component  type  l “ ; UAL$< Ns  ) ;  “  3 " 

12180  LET  Pi x 1 ( 1  .Lent  )-Ns 

12190  NEXT  Lent 

12200  ELSE 

12210  LET  Is-Lxtnt 

12220  FOR  Lent- I  TO  9 

12230  DISP  1st  “layers  remain,  how  many?  or  what  fraction?'’; 

12240  DISP  “  belong  with  type  r “ »UAL$< Lent  ) t " 3 “ ; 

12250  INPUT  “".How 

12260  IF  How<;0  THEN  LET  Lcnt-10 

12270  IF  how/0  THEN 

i  2280  IF  how*,  i  ThEn 

i2290  lET  ns-Inh How»Lxtnt*.5  > 

i 2300  ElSE 

i23i0  lET  Ns-InT(How*.5 > 

i 2320  EnD  IF 
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2330  Fuft  Ftx-Is  TO  Is-Ns  STEP  -i 

2340  IF  F 1  x .> 0  THEN  LET  Pixi<  i .Fix >-Lcnt 

2350  NEXT  Fix 

2360  LET  Is-Is-Ns 

2370  END  IF  i  H oui/0  test 

2360  IF  Is<.i  THEn  Lent* i 0 

2390  NEXT  Lent 

2400  FOR  Lcnt-i  TO  Lxtnt 

24i0  LET  ns-i +IN n RND»Lxint  ) 

2420  IF  LentONs  HnD  Nsx-Lxtnt  THEN 

2430  LET  Swap-r ix i u  .Lent > 

2440  LET  Pixit  i .Lent  )-Pixit  i  ,Ns > 

2450  lET  Fix  it  i  ,ns  >»Swap 

2460  EnG  IF 

2470  NEXT  Lent 

2460  END  IF  1  Strai  test 

2490  IF  Stack  TmEn  i  rest  paraiiei  stacking 

2500  FOR  Lcnt-i  TO  Lxtnt 

2510  FOR  Fix-2  Tu  Lxtnt 

2520  LET  Pixi< Fix  .Lent >-Pixi< i  .Lent  ) 

2530  NEXT  Fix 

2540  NEXT  Lent 

2550  ELSE  i  rest  by  series  stacking 

2560  FOR  Lcnt-Lxtnt  TO  i  STEP  -i 

2570  lET  Suiap-Pixi<  i  .Lent  ) 

2560  FOR  Fix-Lxtnt  TO  i  STEP  - i 

2590  LET  PixHLcnt  .Fix  )-5wap 

2600  NExT  Fix 

26  i0  NEXT  Lent 

2620  END  IF  >  Stack  test 

2630  SuBEND 

2640  i  ii  j[  II  ii  ji  ii  ii  ]i  Ii  ii  Ii  ii  It  ii  II  ii  li  It  It  li  II  II  II  It 
2650  506  Pixi2d_cbox 

2660  i  suproutine  to  make  concentric  Poxes 

2670  COn  /Pass/  Reiay  1  for  sharing  to  suds 

2660  COn  /Pixel/  Cndr*C60i  ,Dhar*i60I  .INTEGER  Lxtnt ,P ix 1 1 i : 1 60 , i : 1 60  / 

2690  INTEGER  Cntr  ,Snei i  .Fix  .Lent  , Peri  .Cntrk  .Past 
2700  REOIri  Pixit  i  ! Lxtnt  ,i  s Lxtnt  ) 

27  i 0  INPUT  “Start  where7  corner-0  center- 1  “  ,Cntr 

2720  IF  Cntr <0  THEN  STOP 

2750  LET  Cntr-BITt Cntr  ,0 > 

2740  IF  Cntr  THEN 

2750  lET  Cntrk»5HIFTt Lxtnt , i  > 

2760  ELSE 

2770  LET  Cntrk-Lxtnt 

2780  END  IF 

2790  LET  Past-i  I  aefauit  assignment 

2800  FOR  Shell-i  TO  Cntrk 

26i0  lET  Fix-0 

2820  LET  Peri -ShIFT t  Snei i  i  J-  i  >  i / 4  per imeter-2»5nei i- t 

2830  DISP  nt  concentric  box  she! i  ; She i 1 »  containing  '  » 

2840  DISP  SHIFT t  Per i  , -SHI FT  tCntr,-i  )M  pixels,  assign; 

2650  INPuT  '  component  type  I i,.. 53  .Fix 

2660  IF  Fix<0  ThEN  STOP 
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12670  IF  Fix-0  THEN 

i 2080  lET  Fix-rasi 

i 2660  ElSE 

i2S00  LET  Fast-Fix 

i26i0  Enu  IF 

i  2 020  PRINT  Shei 1 ; i Snei 1 s  containing  i SHIFT! Peri  .-SHIFT! Cntr  i  J a 

12930  PRINT  pixels  is  assigned  component  type,  E  i: ;  v/fiLSi  Fix  ) ; "  ]  :  i  :  from 

1 2040  IF  Cntr  TmEn  PRInT  ::  center. 

12950  IF  nuT  \Cntr;  THEN  PRINT  :  corner. 

i2560  FOR  Lcnt-i  TO  SniFTi  Snei  i  ,-Cntr  > 

i2970  IF  Cntr  ThEn 

1 2 880  lET  Pixm Cntrk-Snei i+l  .Cntrk-Shei i+Lcnt /-Fix 

i2990  lET  PixitCntrk+Sneii .Cntrk+Sneii-Lcnt+i  )-Fix 

i  3000  LET  P ix i <  Cntrk-Snei  i+Lcnt  .Cntrki-Snei i /-Fix 

i30i0  lET  Pixn  Cntrk+Snei i-Lcnt-t-i  ,Cntrk-Shei  1-m  /-Fix 

i 3020  ELSE 

i 3050  lET  Fix i i Snei i  .Lent /-Fix 

i  5040  lET  rixii Lent  ,5heii/-Fix 

13050  EnD  IF 

i3060  nEaT  Lent 

i3070  nEaT  Sneii 
13060  SuBEnO 

i  3090  •  K  iC  a  H  n  ][  ]i  1 1  3i  ii  H  II  ji  3C  ]i  H  i[  a  ]t  It  Ii  ][  3i 
i3j00  REn  H  sut/routine  to  make  a  Pixel  grid  witn  an  ellipse  inclusion 
i5ii0  REn  lmDeaded  in  a  nost.  The  eiiipse  component  is  always  tne  smaiier 
1 3 1 20  REn  of  tne  two  componnents  and  thus  is  symmetric  witn  respect  to  tne 
1 3 1 30  REn  fractional  volume  filling  factor.  Since  tms  is  limited  to.oniy 
i  3 1 40  REn  two  components  tnen  tne  volume  of  the  2nd  is  used  as  the  variaoie. 
i3i50'Cun  /Pixei/  Cndr*i60]  ,uhdr!i>C60]  .INTEGER  Lxtnt  .Fix  it  I  : 1 80  ,  l : 1 60  / 
i  3 i 60 ' Integer  Hparm  ,noof  ,6de  ,Pof  ,Ptst ,Prq ,Sqrs ,Swup ,Swpxs ,Ap , rp 


a 
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Appendix  II 

Plot  Program  for  Numerical  Solutions 


10  !&&&&&&&&&&&&&&&&& 

20  1  “PL0T_NWK“  is  a  program  designed  to  input  dielectric  data  from 

30  !  network  simulations  and  plot  it.  July  1990 

40  COM  /Memr/  Uctrs( 1 : 505 , 1 : 5  )  ,Hdra*C 80 1  ,Hdrb* I 80 ]  ,Sf i ISC B4 ]  . INTEGER  Dmx  ,Dmy 
50  INTEGER  Kif.Xpt  ,Ypt  .Pit  ,Xlmn,Ylmn,Rpt  ,Sub j .Mach  ,Ef  f  ,Cpy 
60  DIM  ClbSIt]  ,HlbSI40]  ,Xlb$I40] ,Ylb*t40] 

70  LET  Kif-0 

00  CLEAR  SCREEN 

90  PRINT  RPT*<“  "  ,24  )  j  “  A  GRAPH  IS  WORTH" 

100  PRINT  RPTS< "  * , 2 4  ) ; " M  A  N  Y  DATA  POINTS" 

110  PRINT  RPT*< “  “  ,27  )jDATE$<  TIMEDATE  >i *  at  ‘ tTIMESI TIMEDATE  ) 

120  PRINT 

130  !*••>  Data  preparation 

140  INPUT  “  Indicate  data  source:  0-mternal  via  COM  /Memr/  1-file  ",Kif 
150  IF  K i f <0  THEN  STOP 

160  IF  KiP-1  THEN 

170  INPUT  "  Name  the  data  source" “ fi le" *"  ,Sf i 1* 

180  PRINT  “  Data  comes  from s Sf i 1* ; “ “ " : " 

190  ASSIGN  SSourc  TO  Sf ilSjFORMAT  OFF 
200  ENTER  0Sourc jHdra*  .Hdrb*  ,0mx ,Dmy 

210  REDIM  Uctrs< 1 :Dmx  ,  1 :0my > 

220  ENTER  SSourc «Uctrs< * > 

230  ASSIGN  flSourc  TO  * 

240  ELSE 

250  PRINT  "  Data  internal:" 

260  IF  Dmx>0  AND  Dmy>0  THEN  REOIM  Uctrs< 1 :Dmx  ,  1 :Dmy  ) 

270  END  IF 

280  PRINT  HdraS 

290  PRINT  HdrbS 

300  PRINT  “  Data  stored  as  arrayt  1  :“?Dmxi“,  1  :"iDmyi"  )"j 
310  IF  Dmx< 1  OR  Dmy<1  THEN  PRINT  “  ???“ 

320  IF  Dmx>0  AND  Omy>0  THEN  PRINT 
330  LET  Kif-0 

340  INPUT  “  Type:  0)  to  skip  listing  1)  to  list  array", Kif 

350  IF  Kif <0  THEN  STOP 

360  IF  Kif-1  THEN 

370  FOR  Xpt-1  TO  Dmx 

380  FOR  Ypt-1  TO  Dmy 

390  PRINT  Uctrs(Xpt  ,Ypt  )  , 

400  NEXT  Ypt 

410  PRINT 

420  NEXT  Xpt 

430  WAIT  .5 

440  END  IF 
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450  LET  Subj-0 

460  INPUT  "Plot  manner?  defined  iater-0  ser<->pr 1 1-1  ,  diel-2  LOGS-3"  ,Subj 
470  IF  Subj<0  THEN  STOP 

480  IF  Subj-0  THEN  LET  Subj-1  1  Oefaults  to  Ser ies<->Paral lei  plot 
490  LET  Xlmn-0 

500  DISP  "Enter  column  vector  1 . . " i UAL$( Dmy > ; “  for  X  axis  "i 
510  INPUT  “data'?"  .Xlrnn 

520  IF  (Subj-1  OR  Subj-2)  AND  Xlmn-0  THEN  LET  Xlmn-2  !  Default  in  ser<->prll 

530  LET  Aaa-MAXREAL 

540  LET  2 zz-MINREAL 

550  FOR  Xpt-1  TO  Dm.x 

560  LET  Tst-Uctrs( Xpt  ,X lmn  ) 

570  IF  Tat<Aaa  THEN  LET  Aaa-Tst 

580  IF  Tst>Zzz  THEN  LET  Zzz-Tst 

590  NEXT  Xpt 

600  IF  Subj-1  THEN 

610  LET  X loui-0 

620  LET  Xhigh»1 

630  ELSE 

640  DISP  "Minimum  of  X  axis  data  is  " lAaai",  what  shall  be  minimum  for"; 

650  INPUT  "  plot?"  ,Xlow 

660  DISP  "Maximum  of  X  axis  data  is  "tZzzs",  what  shall  be  maximum  for"; 

670  INPUT  "  plot?"  .Xhigh 

680  END  IF 

690  IF  Subj-3  THEN 

700  LET  X low-LGT <  ABS( X  low  > ) 

710  LET  Xhigh-LGT ( ABS<  Xhigh  ) ) 

720  END  IF 

730  LET  Xspan-Xhigh-X low 

740  LET  Xavg-( Xhigh+X low  )• . 5 

750  IF  Subj-3  THEN  PRINT  "  LOG 

760  PRINT  “  X  axis  of  plot  to  range  from  “»Xlow;"to  "iXhigh;"." 

770  LET  Ylmn-0 

780  DISP  "Enter  column  vector  1 . . " ;UAL$( Dmy  "  for  Y  axis 
790  INPUT  “data'?"  ,Ylmn 

800  IF  Subj-1  AND  Ylmn-0  THEN  LET  Ylmn-4  !  Default  in  ser<->prll  case 

801  IF  Subj-2  AND  Ylmn-0  THEN  LET  Ylmn-3  !  Default  to  permittivity  case 

810  LET  Aaa-MAXREAL 

820  LET  Zzz-MINREAL 

830  FOR  Xpt-1  TO  Omx 

840  LET  Tst-Uctrs< Xpt  ,Ylmn  ) 

850  IF  Tst<Aaa  THEN  LET  Aaa-Tst 

860  IF  Tst>Zzz  THEN  LET  Zzz-Tst 

870  NEXT  Xpt 

880  IF  Subj-1  THEN 

890  LET  Ylow— 1 

900  LET  Yhigh-1 

910  ELSE 

920  DISP  “Minimum  of  Y  axis  data  is  “;Aaa;",  what  shall  be  minimum  for"; 

930  INPUT  *  plot?", Ylow 

940  DISP  "Maximum  of  Y  axis  data  is  *iZzz»",  what  shall  be  maximum  for"; 

950  INPUT  *  plot?" ,Yhigh 

960  ENO  IF 

970  IF  Sub i-3  THEN 
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980  LET  Ylow-L6T ( ABS( Ylow  ) ) 

990  LET  Yhigh-LGT< ABS< Yhigh > > 

1800  END  IF 

1010  LET  Yspan-Yhigh-Ylow 
1020  LET  Yavg-< Yhigh+Ylow  )*.S 
1030  LET  Eff-0 

1040  IF  Subj-1  THEN  c.f 

1050  INPUT  "Overlay  an  effective  medium  [EM]  curve  for  2D?  no-0  yea-i  ,trf 

1060  IF  Eff<0  THEN  STOP 

1070  LET  Eff-BIT(Eff.0> 

1080  IF  Eff  THEN  INPUT  "Enter  permittivity  ratio  for  EM  curve  .Kern 
1090  IF  Rem<0  THEN  STOP 

1100  IF  Rem-0  THEN  LET  Eff-0 

1110  ENO  IF 

1120  IF  Subj-3  THEN  PRINT  "  L06  “i  .  . 

1130  PRINT  "  Y  axis  of  plot  to  range  from  “5Yloui"to  ";Yhighi  . 

1140  SELECT  Sub j 

1 1 50  CASE  -0 

1160  INPUT  "Title?"  .HlbS 

1170  IF  HlbS-""  THEN  Hlb*-Hdra*t 1 ,40 3 

1180  INPUT  “X  axis  label?"  ,Xlb* 

1190  INPUT  “Y  axis  label?", YlbS 

1200  CASE  -1 

1210  LET  Hlb*-“RANDOM  2D  NETWORKS* 

1220  LET  XlbS-“Uolume  fraction  of  the  larger" 

1230  LET  Ylb$-"Averaging  Parameter" 

1240  CASE  -2 

1250  LET  Hlb*-”RANDOM  20  NETWORKS" 

1260  LET  Xlb*-"Uolume  fraction  of  the  larger" 

1270  LET  Ylb*-"Permittivity“ 

1 280  CASE  -3 

1290  INPUT  "LOG  Title?", Hlb* 

1300  IF  Hlb$-““  THEN  HlbS-HdraSC 1  ,40 3 

1310  INPUT  “LOG  X  axis  label?', Xlb* 

1320  INPUT  “LOG  Y  axis  label?"  ,Ylb* 

1330  END  SELECT 

1340  INPUT  "  Plot  data  0)  by  points  1>  by  lines', Pit 
1350  IF  Plt<0  THEN  STOP 
1360  IF  Plt-0  THEN 

1370  INPUT  "  Type  a  single  letter  character  for  the  plot  points  .Lib* 

1380  INPUT  "  Character  size  in  I  of  Xs  of  the  graph  width?  le  1,2,..  ,Csz 

1390  IF  Csz-0  THEN  LET  Csz-1  !  Defaualt 

1400  END  IF 

1410  IF  Mach-0  THEN  »  m 

1420  INPUT  "Hardcopy  of  plot?  0-none  1-to  printer  2-to  plotter  .Placrt 

1430  IF  Mach<0  THEN  STOP 

1440  END  IF 

1450  IF  Mach-2  THEN  „ 

1460  INPUT  "Max  speed  of  plotter  pen?  (*1-20,  in  cm/s)'  .Speed 

1470  IF  Speed<0  THEN  Stop 

1480  IF  Speed-0  THEN  LET  Speed-10 

1490  END  IF 

1500  WAIT  1 

1510  CLEAR  SCREEN 
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1520 

KEY  LABELS  OFF 

1530 

!**#>  Begin  plotting 

1540 

GINIT 

1550 

PLOTTER  IS  CRT  INTERNAL" 

1560 

GRAPHICS  ON 

1570 

FOR  Cpy-0  TO  (Mach>0) 

NEAR  FULL  SIZE  SCREEN 

1580 

i VIEWPORT  20.120,15,85 

t 

1590 

IF  Sub j “3  THEN 

1600 

VIEWPORT  40,95,35,90 

1610 

ELSE 

1620 

VIEWPORT  30,1 10,35,90 

• 

1630 

END  IF 

1640 

WINDOW  Xlow  .Xhigh ,Ylow .Yhigh 

1650 

IF  Plt-0  THEN 

• 

1660 

FOR  Xpt-1  TO  Qmx 

1670 

LET  Xx-Vctra( Xpt  ,Xlmn ) 

1680 

•IF  Subj-1  THEN  Xx-1-Xx 

I 

IF  COMPLEMENTING  NEEDED 

1690 

LET  Yy-Vctrs( Xpt ,Ylmn> 

1700 

IF  Subj-3  THEN 

1710 

LET  Xx-LGT(AB$<Xx ) ) 

| 

1720 

LET  Yy-LGT< A8S< Yy  )  ) 

1730 

END  IF 

1740 

MOVE  Xx ,Yy 

1750 

CS1ZE  C3Z..5 

1760 

LORG  5 

1770 

IF  Clb$-"“  THEN 

1780 

LET  ClbS-“+“ 

1790 

IF  SubjOl  ANO  Onx<  1 00 

THEN  LET 

Clb*-VAL$( Xpt  ) 

1800 

END  IF 

1810 

LABEL  Clb* 

1820 

NEXT  Xpt 

1830 

ELSE 

1840 

!***>  curve  higher  by  a  std 

deviation 

1850 

i  FOR  Xpt -2  TO  Onx-I 

1860 

i  LET  Xx-Vctra(  Xpt  ,1  > 

1870 

i  LET  Yy-Vctra(Xpt  ,3 >+Vctra< Xpt  ,6 > 

1880 

i  LINE  TYPE  3 

1890 

i  PLOT  Xx  , Yy 

1900 

'  NEXT  Xpt 

1910 

•  PENUP 

1920 

•***>  central  curve 

% 

1930 

FOR  Xpt-1  TO  Dmx 

1940 

LET  Xx-Vctrs( Xpt  ,X Inn  ) 

i  IF  COMPLEMENTING  NEEDED 

1950 

i IF  Subj-1  THEN  Xx-1-Xx 

I960 

LET  Yy-Vctrs( Xpt ,Ylnn ) 

1970 

IF  Subj-3  THEN 

1980 

LET  Xx-LGT <  ABS<  Xx  ) ) 

1990 

LET  Yy-LGT ( ABS<  Yy  ) ) 

2000 

ENO  IF 

2010 

LINE  TYPE  1 

2020 

PLOT  Xx.Yy 

2030 

NEXT  Xpt 

2040 

PENUP 

2050 

END  IF 
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!  overlay  effective  medium  curve 


2060  IF  Eff  THEN 

2070  FOR  Xpt-1  TO  200 

2080  LET  Xx-Xpt/201 

2090  Aaa-(Rem-) )*<2*Xx-1  ) 

2100  LET  Oyf-.S*< Aaa+SQR< Aaa*Aaa+4*Rem  )  > 

2110  LET  Yy-FNWnr(Rem,Xx  .Oyf  ,Zzz> 

2)20  ip  ABS( Zzz  )> .000001  THEN  PRINT  "  Warning!  Ser'«>Prll  difficulties 

2130  LINE  TYPE  1 

2140  PLOT  Xx ,Yy 

2150  NEXT  Xpt 

2160  END  IF 

2170  IF  Subj-1  THEN 

2180  AXES  .05, .1  ,0,-'  .2,1 

2190  AXES  .05..  1  ,1  ,1  .2,1 

2200  ELSE 

2210  FRAME 

2220  END  IF 

2230  !***>  set  labels 

2240  CLIP  OFF 

2250  L0R6  4 

2260  MOVE  Xavg ,YhigH+.08*Yspan 

2270  CSIZE  4,. 6 

2280  LABEL  Hlb$ 

2290  LORG  4 

2300  MODE  Xavg  ,Ylow- . 1 6*Yspan 

2310  CSIZE  4, .6 

2320  LABEL  Xlb$ 

2330  MODE  Xlow-( . 1+(Subj-1  >* . 04  )*Xspan , Yavg 

2340  DEG 

2350  LOIR  90 

2360  LORG  4 

2370  LABEL  Ylb$ 

2380  IF  Subj-1  THEN 
2390  CSIZE  2. 4,. 6 

2400  LORG  1 

2410  MODE  X low-. 092 *X span ,Ylow 

2420  LABEL  "Series-1 ike“ 

2430  LORG  7 

2440  MODE  X low-. 092 *X span  ,Yhigh 

2450  LABEL  "Parallel-like" 

2460  LDIR  0 

2470  LORG  5 

2480  CSIZE  2. 2, .6 

2490  LET  Zzz-Ylow- . 04 *Y span 

2500  FOR  Xgf-Xlow  TO  Xhigh  STEP  .1 

2510  MODE  Xgf- .01 ,Zzz 

2520  LABEL  USING  "O.D"iXgf 

2530  NEXT  Xgf 

2540  LET  Aaa-X low- .05*X span 

2550  FOR  Ygf-Ylow  TO  Yhigh  STEP  .2 

2560  MODE  Aaa,Ygf+.01 

2570  ILDIR  90 

2580  LABEL  USIN6  “DO.D" »Ygf 

2590  NEXT  Ygf 
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2600  END  IF 

2610  LORS  1 

2620  LDIR  0 

2630  MOVE  X low- . 1 9»Xspan  , Y low- . 30*Yspan 
2640  CSIZE  2, .6 

2650  LABEL  HdraS 

2660  MOVE  Xloui-.  1  9»Xspan  ,Ylow-.35*Yspan 

2670  LABEL  Hdrb$ 

2680  SELECT  Mach 

2690  CASE  -1 

2700  DUMP  GRAPHICS  CRT 

2710  WAIT  2 

2720  LET  Cpy-2 

2730  CASE  -2 

2740  GINIT 

2750  PLOTTER  IS  705,"HPGL“ 

2760  GSENO  ''VS''&VAL$<  Speed  ) 

2770  END  SELECT 

2780  NEXT  Cpy 
2790  KEY  LABELS  ON 
2800  PEN  0 
2810  END 

2820  DEF  FNWnr( Drat io  ,Frpx  ,Din  .Pcterr > 

2830  !  *  *  *  >  Object  of  function  to  find  alf  satisfying  1 “SUM  of ( vol *diel Aal f  ) 

2840  i**#>  where  vol-fract ional  volumes 

2850  !***>  diel-( species  permittivity  )/( composite  permittivity) 

2860  !*•*>  alf-  constant  exponential  between  -1  &  +1  ( ser ies<->paral lei  ) 

2870  INTEGER  Spk .Arnd 

2880  IF  Din-0  OR  Frpx>-1  THEN 

2890  LET  Alf 1-0 

2900  LET  Pcterr-1 .E-99 

2910  PRINT  "  ..  series  <->  parallel  factor  at  or  beyond  limits'' 

2920  ELSE 

2930  IF  Dratio>0  AND  Frpx>0  ANO  Frpx<1  THEN 

2940  LET  Gratio-L0G< 1 /Din  ) 

2950  LET  Gsum-( 1 -Frpx  )*Grat io 

2960  LET  Gdev-( 1 -Frpx  )*Grat io*6rat io 

2970  LET  Grat io-L0G< Drat io/Din  ) 

2980  LET  Gsum-6sum+Frpx*Grat io 

2990  LET  Gdev-Gdev+Frpx*Grat io*Grat io 

3000  END  IF 

3010  LET  A1 f 1 --2»Gsum/Gdev 

3020  LET  6sum-Gsum+Gdev  !  Initial  guess  of  sum 

3030  LET  Try-1  I  Optimize  sign  in  iterations 

3040  LET  A1 f 0-0 

3050  LET  Arnd-1 

3060  WHILE  ABSIAlf 1-Alf0)>. 0000000001 

3070  LET  Gsum-0 

3080  LET  Gdev-0 

3090  LET  Gratio-1/Din  !  relative  permittivity  to  composite 

3100  LET  Gwk-0  i  a  term  0  species  to  sum 

3110  IF  Grat io< >0  AND  Frpx>0  AND  Frpx<1  THEN 

3120  LET  Gwk-( 1 -Frpx  )»6rat io" A1 f 1 

3130  LET  Gsum-Gwk  !  sum  function  (as  described) 


45 


!  1st  derivative 


3140 
3150 
3  <60 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3370 
3380 
3390 


LET  bdev*LUb( bratio >*6wk 
END  IF 

LET  Grat io-Drat io/Din  !  relative  permittivity  to  composite 

LET  Gwk-0  !  a  term  9  species  to  sum 

IF  Grat io<>0  ANO  Frpx>0  AND  Frpx<1  THEN 
LET  Guik-Frpx*Grat io” Ai f  1 

LET  Gsum-Gsum+Gwk  !  sum  function  (as  described) 

LET  Gdev“Gdev+LOG< Grat io  )*Guk  !  1st  derivative 
END  IF 

LET  Pcterr-Gsum 

LET  Gsum-Al f 1 *( Gsum-1  >/Gdev 

LET  Gdev-Al f 1 *A1  f 1 -Try*6sum 

IF  Gdev<0  THEN  Gdev-Gdev+2*Try*Gsum 

LET  Alf2-SGN< Alf 1  )*SQR(Gdev> 

IF  Arnd  MOD  2-0  AND  A8S< Al f 2-Ai f 1  >>ABS< Al f 1 -Al f0 >  THEN  Try— Try 
LET  Alf0-Alf1  !  bumping  iterations 

LET  Al f 1 -Al f 2 

LET  Arnd-Arnd+1  !  incrementor 

END  WHILE 

IF  Al f 1 < >0  AND  Pcterr<>0  THEN 
LET  Pc  t err- 1 00* (  ABS<  Pc  t err  )*  (  1  /  A‘l  f  1  >-1  ) 

END  IF 
END  IF 
RETURN  Al f 1 
FNEND 

!  H  ][  It  K  K  H  K  H  H  It  It  It  It  It  It  It  It  It  It  It  It  It  It  It 
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MANDATORY  DISTRIBUTION  LIST 
CONTRACT  OR  IN-HOUSE  TECHNICAL  REPORTS 


101  Oefense  Technical  Information  Center* 

ATTN:  DTIC-FDAC 

Cameron  Station  (Bldg  5)  (*wote:  Two  copies  for  otic  will 

Alexandria,  VA  22304-6145  *  *««*  ^o«  stinfo  Office.) 

483  Director 

US  Army  Material  Systems  Analysis  Actv 
ATTN:  DRXSY-MP 

001  Aberdeen  Proving  Ground,  MO  21005 

563  Commander,  AMC 
ATTN:  AMCDE-SC 
5001  Eisenhower  Ave. 

001  Alexandria,  VA  22333-0001 

609  Commander,  LAB COM 

ATTN:  AMSLC-CG,  CD,  CS  (In  turn) 

2800  Powder  Mill  Road 
001  Adel  phi,  Md  20783-1145 

612  Commander,  LABCOM 
ATTN:  AMSLC-CT 
2800  Powder  Mill  Road 
001  Adel  phi ,  MD  20783-1145 

680  Commander, 

US  Army  Laboratory  Command 
Fort  Monmouth,  NJ  07703-5000 

1  -  SLCET-DD 

2  -  SLCET-DT  (M.  Howard) 

1  -  SLCET-DB 

35  -  Originating  Office 

681  Commander,  CECOM 

R&D  Technical  Library 

Fort  Monmouth,  NJ  07703-5000 

1  -  ASQNC-ELC-I-T  (Tech  Library) 

3  -  ASQNC-ELC-I-T  (STINFO) 

705  Advisory  Group  on  Electron  Devices 
201  Varick  Street,  9th  Floor 
002  New  York,  NY  10014-4877 
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205 

Oi rector 

603 

Cdr,  Atmospheric  Sciences  Lab 

Naval  Research  Laboratory 

LAB COM 

ATTN:  CODE  2627 

ATTN:  SLCAS-SY-S 

001 

Washington,  OC  20375-5000 

001 

White  Sands  Missile  Range,  NM  88002 

221 

Cdr,  PM  JTFUSION 

607 

Cdr,  Harry  Diamond  Laboratories 

ATTN:  JTF 

ATTN:  SLCHO-CO,  TO  (In  turn) 

1500  Planning  Research  Orive 

2800  Powder  Mill  Road 

001 

McLean,  VA  22102 

001 

Adel  phi ,  MD  20783-1145 

301  Rome  Air  Development  Center 

ATTN:  Documents  Library  ( TILO) 
001  Griffiss  AFB,  NY  13441 


437  Oeputy  for  Science  &  Technology 
Office,  Asst  Sec  Army  ( R&O) 

001  Washington,  DC  20310 

438  HQOA  (OAMA-ARZ-O/Dr.  F.O.  Verderame) 

001  Washington,  DC  20310 

520  Oir,  Electronic  Warfare/Reconnaissance 

Surveillance  and  Target  Acquisition  Ctr 
ATTN:  AMSEL-EW-D 

001  Fort  Monmouth,  NJ  07703-5000 

523  Oir,  Reconnaissance  Surveillance  and 
Target  Acquisition  Systems  Directorate 
ATTN:  AMSEL-EW-DR 

001  Fort  Monmouth,  NJ  07703-5000 

524  Cdr,  Marine  Corps  Liaison  Office 

ATTN:  AMSEL-LN-MC 

001  Fort  Monmouth,  NJ  07703-5000 

564  Oir,  US  Army  Signals  Warfare  Ctr 

ATTN:  AMSEL-SW-OS 
Vint  Hill  Farms  Station 
001  Warrenton,  VA  22186-5100 

602  Oir,  Night  Vision  &  Electro-Optics  Ctr 

CECGM 

ATTN:  AMSEL-NV-0 

001  Fort  Belvoir,  VA  22060-5677 
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